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ABSTRACT 

The pa>es:ent, voi>k Is eri Co drvelope e PC b^sed 

irxCex'ecdve peoke^e for> px»edicdri^ dispex'siori olT pollutCenCs 

emiCCed fiorn potnC soux’CeCs!) in Che eCnxosphene The peickie^e is 

wniCCen in TURBO PASCAL vex^sion 4 0 Co pnovide ^nephlc visiiad arid 
Co aohdeve user- macMne inCenacCion The selecCion is achieved 
Chnou^h a scneen meriu usin^ annow keys Che keyboa 3 >d 

The package also provides faoiliCies Co have access Co any ’wond 
processing sofCwax^e CVORDSTAR 2000 cunnenCly implemented> Co 
edlC/^cr-eaCe inpuCy'ouCpuC files aC appropriaCe levels 

The menus ai^e divided inCo Cvo levels and Che menu 

applicable aC a parCiculax sCa^e of operaCion is displayed The 

selecCioris made ax^e displayed simulCaneously indicaCln^ Che paCh 
Chrou^h which an individual Cenminal selecCion is idendfied When 
a selecCion is being- made OPTION ON is displayed while aC Che end 

of selection READY mmy Is displayed which per^mlCs conClnuirig with 
Che cux’renC selection^ making addlClonal selection^ printing Che 
screen display ox Cei^mlnaClon of Che package as may be desired 
The package inlClallzaCion and execution starts by giving command 
INDISPOL On line HELP and GUIDE are available for package 

operation and making appxopriaCe selections x^especdvely 

The simulation models employed for dispersion of 

pollutants ax^e classified Into five categories^ UNIS for single^ 

steady and continuous point source; MULTYS fox' steady and multiple 
continuous point sources; PUFFS foX' fluctuating point emissions in 
calmwind condition; SEGMNT for fluctuating point emissloxis in 

transport condition^ and PUFSEG for fluctuating point emissioris In 
llghtwlxid conditions UNIS and MULTYS assume statlonairy and 
homogeneous atmosphere and use Gaussian approach PUFFS^ SEGMNT 
and PUFSEG consldex' non-- stationary and inhomogeneous atmosphere 
and use extended Gaussian appx-oach 

The package structure provides facilities for both 

digital as well as graphical output Howevex ^ the current version 

suppox'ts gx aphlcal facilities for UNIS only 


KEYWORDS 

Alt' pollution^ atmospheric dispex^slon> Gaussian 

diffusion^ non-stationar’y atmosphere^ inhomogeneous atmosphez-e^ 
puff^ segmented plume^ mixed puff-segment approach^ shox't-term 
emlssiorxsj, unsteady-state conditions^ computer software^;, MS DOS^ 
Turbo pascal^ intex active package^ interactive graphics^ personal 
computer's^ IBM compatible 



1 INTRODUCTION 


An eff'ec-tive pnot.ect,ion of oux- envix^onmont is highly 

dep€?rid^riL on ppoclse inf onm^i'Lion on i-ts s'LaL'to aarid developmerit. In 
Ihiis way> Gomput,<&x» use In envix^onmenfel agencies hes been ^ednln^ 
mone arid more imporfance However ^ as bl sclentrlflo applicatrion 
fields envlronmenbad compubin^ was only of Umibed research 
Inberesb bill now Just recently an Incx'easlng^ research Interest 

cari be observed as evident by intensified national and 
internatinal publishing Ce ^ ENVIROSOFT joux*nal> and recent 

conferences Including the on^oin^ ENVIROSOFT-serles 

In general ^ different functional application fields for 
computers in environmental management can foe stated as follows 

Data CO llectioriy measuring <e ^ measuring of air 

pollutants in air quality monitoring systems > 

Data storage and retrieval Ce ^ envirorimental research 

project or literature databases with fle>dfole search furictloris> 

Data analysis <e ^ estimation of emissioris from 

computerised emission registers J 

Modelling arid simulation Ce ^ propagation calcirlations 
for air or water pollutantsJ 

Grap flics Ce ^ land use mapping ^ 

A survey report CPa^e^ 1988> Indicate that computers are 
mainly employed fox data processings C77% applications) while 
modelling and simulation applications are only 12% On the system 
tecfirdcal slde^ the study CPage^ 1988> revealed a dominance of 
main frame programs Cmore than 50%) which are run on Inhouse 
machilnes or in community centres The PC-use in environmental 
computing^ however^ is catching up Cmore than 2D%> The 
predominant programming lanauge is still FORTRAN C45%)^ modern 
lanauges such as PASCAL or the database languages ar^-e not yet very 
common C10%) The statistics of the number of users per programme 
shows that environmental softwar^e is only used by minority of the 



in envlponmen'tal orrices <145% of -tho pno^nammee: wi“th maximum 
5^ only S% wi-th mor«e than 10 r-e^uiax* us:ens:> This may bo duo to 
the availability at limited number- of PC based user friendly or 
interactive packa^es^ particularly in the field of modelling and 
simulation 


Modelling and simulation of envix onmental pxoblems hold 
a lon^ tradition <lFienklei arid Goodall, 1978> The software 
techniques used in envir'onmental modelling ^ however^ are not 

always on a satisf actor-y level fx^om the usei point ol view User 
friendly or interactive software cari offer powerful concepts for 
supporting desig^n^ implementation and operation of models which 
axe not yet suf f Icient ly used in envix onmental modeiiin^ These 
&of twax'e concepts provide a hig^h decree in flexibility and use 
adapatlon in xunnin^ envix orimental analysis methods and models 

Intex active modelling and simulation systems ai^e mainly concerned 

with computer^ support of the entire process of model building and 

simulation^ fx-om piarinin^^ to buildln^^ and documentation of 
models up to ^rapWcal presentation Thx’ou^h a hi^h decree of 
interactivity ori all system levels interactive C^raphical> 

modelling dixectly at the display^ interactive control of 
simulation flow as well as interactive preparation of simulation 

r-esults is offered The aim is a M^h decree of friendliness in 

computer modelling partlcuiaxdy fox* usex'S with little ox* no 
computex’ expexience The main asset of this modelling tool is its 

object oriented user interface taking advantage of modern 

graphical featui*es such as. windows^ pull down menus^ etc 

Concepts ot intex'active modelling and simulation aiT'e not 
only lelevant for scientific dlspexslon modelling;, but also tor 
Decision Support Systems in planning and policy making Hence ^ the 
majox thxust in the present study is laid on developing an 
Intei'actlve package fox simulation of dlsper*siori of pollutants in 
the atmosphex'e 



2 STATE OF THE ART ON DISPERSION MODELLING 


2.1 

In envir^orimenlfid aclenoe^ £is In moat, or the flelda of" 
applied acinces^ the mein Inter-eat Ilea in ceuae end elTeot 
r elatlonatilpa oi in oeauel e>^plarietloris: of eventa The desine for* 
ttiia knowledge stems fnom a need to know what actions <Ga\Jises> 
must be taken to achieve a panticulan change Cettect^ in 

envinonmental qxjiality In phllosphen Kanl Poppen^s wonda CPoppei*^ 
1950> "castjol €?ycplanatiiDn.s o/ any euent ne<j>uir‘e tuAO In^r-adti.&nts: 
uniuer-sl statement* which ore hypotheses af the chcrr-aeten af 
natunal lows, and singular- stcrtements^ which per-tcrin to the 
spec^/ic event ^n <gfuestion It ^s from tHiS- set o/ these two types 
o/ stotements thet sirij^ulor statements or spec^/^e predf^ctions 
about events are dedfucec# *' 

It is obvious that atmoaphenic motion plays a major- role 
in air pollution phenomena As such attempts must be made Cl> to 

Tormulate new /evaluate existing unlveraal mathematical state- 

mentCsJ) about the char acter of this role^ and C2> to know how 
these statements in conjuriction with meteorological observations 
ultimately affect both the form of the dispersion model and the 
types of predlctioris that can be made The meanln^fxjd use of air 
pollution models of universal statements as a tool In decision 
making process requires both an undex’standlng of the inherent 

limits of model capabilities and some measure of uncertainty 

associated with a particular model's predictions Thus in 
f or mulatln^/^ev aluatin^ /'selecting a dispersion model consideration 
must be ^Iven to the analysis of model's capabilities^ 

uncertainties and simplicity 

2.2 Basic Framework: for Dispersion Modelllr\g 

It is well known that the dlspex^sion of g^aseous and 
partlciriate material released into the atmosphere is dependent on 
turbulent and convective motions px^esent within the boundax'y 



Layer* The trUr-bulen'L moblori £t]>l£:es: becaubsie of t#he physlcel 

pr*oper*"ties: at t,he maber-lel;, 'the vex^-ticad t.empex*a‘tu3>e st,r*uc‘Lux*e arad 
the effect of ground i-ou^hine&e: elements: whene the dr*a^ on the 
sunface is b a l an ced by momentum transmitted through the layer by 
shearing stresses 

The fundamental problem of turbulent motion Cdiffusion>^ 
both practically and theoretically^ Is to predict the vairlations 
with space arid time of the statistical properties associated with 
the distribution of concentration of contaminant in a cloud 

coritalniri^ a flrilte cfuaritlty of material Althoug^h there have been 
notable success on the theoretical slde^ there ai*e formidable 
difficulties preventing a complete solution So prediction in 
practical cases has stilly unforturiately^ to be made usln^ 

empirical assumptions.^ not Justified by <ox even consistent v/lth!> 

theory TMs is ultimately unsatisfactory Even though the use of 
such assur\ptlons leads to reasonable predictions in some^ even 
many cases ^ the basic physics is avoided Thus there is no way of 
assessing the magnitude of errors involved^ which expexdments show 
can be very lar^e Nor does such a procedure help in establishing 
a sure framework fox* practical predictions^ which is definitely 
the ultimate ^oal of all research in atmospheric diffusion 

The understanding and modelling of plume dispersion 

downwind of point and line sources in turbulent boundary layers 
has been hampex’ed by a lack of e>cperlmental lnfor*mation concex*nln^ 
the terms occuring in the exact trarisport equations for the mean 
concentration^ the scalar fluxes and the scalar* fluctuations 

Development of the subject by way of experiments conducted in the 
atmosphere has been limited by the complexities of atmospheric 
boundary layer flows and;, in addition^ the instrumentation arid 

resources necessary to undertake sufficiently detailed experiments 
have generally not been available 

The problem of turbulent diffusion in the atmosphere has 
not yet been uniquely formulated in the sense that a single basic 
physical model capable of explaining all the significant aspects 

of the problem has not yet been proposed There are three 



and 


al'tar-na'Llve apppoachas based on s'tabis'Lioal^ similai^it^y 
^x-adlenb br-ansl^ei- Iheonles CNema^ 19S7; Robins and Paclcrell^ 1079;, 
Svadas^ 1974} None of these cari be cate^ordcaHy eliminated Tnom 
Gonsldenatlon since each one has Its own importance 

A mathematical model of atmospheric dispersion must 
attempt to simulate the ^ross behaviour of plumes emitted from 
g^xourid ox elevated sources Although several basic aforementioned 
approaches ax^e possible^ usually a number of simplifying 
assumptions ax-e necessary in any case to obtain a mathematically 
tractable solution As a r^esult^ ail these theories tend to the 
same distribution function foi' the pollutarit concentration 
Mathematically this mearis that the concentr atlon in the ci'Osswind 
and vextical directions may be Viox mally^ distributed Hence^ the 
most widely accepted model at pr*esent is the Gaussian plume model 
which does exhibit "normar' dlstxdbutlon as su^^ested by the three 
basic tlieoxles 

The Gausslari diffusion model has achieved considerable 
populaLxdty amon^ people attempting to descxdbe the r»ole of 
atmospheric dispersion The ease of application of this model for 
consex vatlve poUutarits out-wel^hts its shortcomings as a complete 
mathematical description of the diffusion px'ocess Most measxired 
dlspex'slon data fit the Gaussian model reasonably well CBowne^ 
1974^ and published dispersion coefficients have been available 
fox corisldex able length of time CPasqulli^ 1974> 

2 3 Recent Developments In Dispersion Modelling 

Extensive wor*k has been done in recent years on air 
pollution dispersion modelling and as a X'esult large number- of 
publlcatloxis ar-e available However-^ none of these publications 
suggest a f undamentally new appx-oach to the problem of dispersion 
modelling and most of them stem from the basic three appi-oaches 
mentioned in the px-evious section The major contribution is in 
terms of Improved idealization of source^ atmosphereXterraln 
cliar acterlstlcs ^ l<Cx-ogstad and Pattex-sori^ 19SS; 21annetti^ 

1984^ Shelh^ 1978> ox- pr-ovldlng algoritluns to solve dispersion 
equations employing numex-lcal techniques Ce g Enger-^ 1986; van 



S'tijn arid Nieuvst^adt.^ 1986 > or’ developing arialyt^icai s.olu'Lioris of^ 
dil'fuisiori equations; under' e-peclallzed condlLlorifi: <e ^ Shukia arid 
Ghiautiari;^ 19863 In apl'te of t,he afor'ement^ioned improvement's: arid 
exparis;lorifi^ it. haue tieen found fliat, offen t,Ke more complex 
methodologies: pos:s:es:s. only a theoretical Cor poteritial> capability 
of better representing the complexities of the real world CLamb^ 
1984> Zarmetti C1986> after reviewing recent important model 
validation studies concludes tliat 


*’wKen the mocfels or-e opplred ^n an operational^ ’hands o//* manner 
shart-t^rm mcjd&lUn^ £:imulat\ans, ar-er £:ubs:tantiall\^ inacc:urat& with 
errors of a factor- of two in more than 5*0% of tho cases, the 
more comployc mode^llin^ appr-oaahos: do not pr-oxjido a ^ubsttantial 
\mpr-o\j&mcfnt %.n n€fpnod tiding r-oalxty^ oompar-od urith the more sttYnplo 
ones " 


Thus there is a need of additional air pollution 

modelling efforts for impi ovln^ the present simulation 
capabilities and allowing the models to reach the level of 

per f or marice that is expected from them,, especially for regulatory 
appllcatloris since air pollution dispersion models are the only 
tool for inferring a quantitative deterministic relation between 

antliJ'Opo^eriic pollutant emlssoris and ambierit concenti-ations 
Zarmetti C19863 presented a new methodology to simulate complex 
dispersion conditions in both transport arid calmwlnd situations 
while maintaining the simplicity of the basic Gaussian eqxiatlon 
The method is claimed to be computationally cost effective and 
allows a non- stationary^ inhomogeneous representation of 
atmospheric phenomena such as transport, turbulent diffusion^ dry 
and wet deposition^ and first order reaction chemistry This new 
approach^ called "mixed Puff-Segment approach*’^ appears to provide 
an Improved slmulaion tool for practical applications in both 

short-rarige and long-range air pollution dispersion studies,. in 
either flat or complex terrain 

In general, atmospheric dispersion is but one of the 

lax'ge class of phenomena which are composed of a deterministic and 
random element The deterministic component may be modelled with 

ail the precision allowed by experimental input while the random, 
stochastic part is unpr edictable Raridomi“iess in atmospheric 
modelling appear'S in two forms 1 the input error, as the initial 



statue can never' be defined wlbhoub er*r* 03 > and wibh Calmest,> 
ini'iril'te def^lnl*tiori^ 2 bbe t'andom cholcea available bo bhe aysbem 
dirrin^ ibs evolubion 

In any sysbem bher-e exlsbs an unequivocal chain bebveen 
an effecb and bhe cause In many cases ib may be a sbx al^hb 
forward process bo describe bhe link bebween cairse arid effecb by a 
mabhemabical evpnession C m models The fonmula may somebimes be a 
very complex one^ bub bhere is always ari equabion Con a sysbem at 
many equaborisO nelabin^ bo bhe effecb bo bhe cause The chain 
bebween cause and effcb is uninodal In abmosphenlc dispension 
bhis chain is mirlbi- nodal The mox-e sophisbicabed bhe dlsper>sion 

calcuiabion^ bhe mone nodes ane pnesenb A "simple** Gaussian 
xeal-bime plume calculabion conbains ableasb bhe boUowln^ nodes 

CBenar'ie^ i9S7 y 1 bo 4^ sounce specibicabions Cemission 
inbensiby^ concenbnabion^ exhausb veloclby^ bempenabune>^ 5 bo 12^ 
mebeonolo^ical panameber-s Cposiblon cooT'dlnabes^ wind velocibies 

and bheir flucbuabions in bhmee coondinabe dlnecbions^ 
bempenabune^ sbablllby>, 13 arid 14^ bhe dlspension pai'amebens C<y^ 
and pollubanb loss <deposiblon> oi' fonmabion paxameben 

Assuming veny opbimisbically ennon of 3% at each node^ bhe end 
result will be ennoneous by about 12 bo 50% Bub if bhe ennor ab 

each node is a still modest 10%> bhe outcome will be uncenbain by 
about 40 bo 400% CBenanie^ 1987^ And all bhis with a physical 
formula wliich is corisider ed as correct 

In Spite of the af onemenbloned inhenenb Umlbablons one 
has bo resonb bo dispension modelling because bhis is bhe only 
bool for* infennin^ a quanbibabive debenmiriisblc nelabion bebween 
pollubanb emissioris arid ambient concentnabions Thus bhene is a 
need bo continue ain pollution modelling ef fonts for* impnovin^ bhe 
pnesenb simulation capabllibies 



3 SCOPE OF THE PRESENT WORK 


The main ■LHr'Ust, in t,he pr'esenh wox*k htas: been on 
developing an Inbeneobive package which caii be lmpiement,ed[ on any 
IBM PC OT' lbs compablble bo &imvilabe poliubanbr dispex'sal l^nom 
single or mulbiple^ conbinuoufi: or* shonb benm^ ebeady or* 

f lucbuabin^^ sui'lTaoe on elevabed s:ounGe<ls> in sbabionairy on 

non-sbablonany^ homogeneous or* inhomogeneous abmosphene The 
package is aimed ab pnovidlng facllibles bon Cl> ble>dhle sysbem 
of selecbion of slmulablon pr ocedune bhnough which bhe usen can 
nevlgabe wlbh ease^ C23 accepbance of inpub daba bhnough benmlnal 
or a file creabed in bhe users area^ C3> checking of Inpub daba 
and ediblng of inpub flles^, C4y guidelines for selecbion of 

simulablorn models Implemenbed^ <5^ builb-ln procedures for 
sbablllby classif Icablon and evaluablon of dispersion paramebers 
in bhe absence of user specif Icabions^ and display arid sborage 

of oubpub in graphical as well as dlgibal form The developmenb of 
a package which will have af oremenbioned facllibies Involves 

following sbeps 

1 Selecbion of appropriabe slmulablon procedureCs> for 

abmosphei Ic dispersion of pollubanbs 

2 Selecbion of programming language which can supporb required 

gr aphd c capabl libies 

3 Selecbion of proper sofbware envixonmenb and supporbing 

sofbware packages bo facllibabe efflcienb ubllizablon of a PC 

4 Developing package programme algor Ibhm and lbs 

implemenbabion 

Selecbion of appropriabe slmulablon bechndque is a very 
difflculb bask Developmenb in abmospherlc dispersion modelling 
becliniques in bhe lasb bwo decades has been quibe remar'kable Wibh 
the parallel growbh in compubablonal capabilibles^ ib has been 
possible bo define and Implemenb exbremely advanced slmulablon 
techniques In general all Improvemenbs and expansions In air 
pollublon dispersion modelling sbem f X'om bhe bhree basic theories 
(sbabisbicaf similar iby and gradlenb bx'ansfer> ox' bhe Gaussian 



appr»oach CChApter* 2> In splt,^ of so much pr-ogness^ none of Lhe 
exls-tin^ modelUn^/'&imuiet.iori ijecFurilqiies yield neftsonehle 

pnedlcllons: Since edn pollutan't dispex'-slon models ^a-e 'the ordy 

l^ool for infer xdn^ e quent^lt^e-tive ne^ixle-tony epplice-tlons^ 

continuous ef fonts ene made ion ImpT’Ovlng: the pnesent simuietion 
capabilities which would adlow the models to neach the expected 
level of penfonmaLnce These efforts aim at C±> the development and 
application ol more complex and sophisticated methodologies;, 

^ener ally requlrdn^ mox'e advanced meteorological infor^mation^ and 
C2y the impx’ovemerit of the simulation capabilities of x^elatively 
simple current techniques^ mainly using the available 
meteorological information Simulation techniques adopted in the 
cux*rent package employ simulation capabilities of relatively 

simple recently developed techniques (INema and Tare^ 1980; Nema^ 
1987 > Zannettl^ 1986 > The methodology appears to simulate complex 
dispersion conditions in both tx-arisport and calmwlnd situation 
while maintaining the simplicity of the basic Gaussian approach 

Selection of px’-ogrammlng language and software 

envlr onment is a relatively easy task since not many choices are 
available TURBO PASCAL vex'slon 4 0 wldch suppor*ts good graphic 
facilities is utilized in the present package The advantage is 
th^t the programmes cauri be executed fr-om outside the TURBO package 
and hence the pi ovlsloris of MS DOS environment are fully utilized 



4- STRUCTURE OF INDISPOL 


4 1 Gerier'fid P£&c:kJal^e^ Pr^o^r'^amme Logic 

The pr-esenL peckeige fists: been developed keeping In view 
Lhe genenai neiLux'e of Lhe packaige arid e>qpeoLed losage by Lhe end 
usen Since Lhe ulLlmaLe objective of Lhe package is bo pr-ovide 
bhe usen a faclUby bo simuilabe pollubanb concenbX’abion pnoflle 
under var^ylng dispersion condibioris using diffenenb simulabion 
bechmlques ajod bo produce oubpub in digibal and/or gi aphical form 
depending upon one^s JudgemenbXrequlremenb^ bhe algoribhm has been 
developed as a mulbl- branched seiecbion process which logically 
ends Inbo bhe execubion of bhe px^ogramme fox’ selecbed opbion and 
display /sborage of oubpub Ib is also assumed bhab bhe user may 
liave libble ox' no pr ogX'amming knowledge Thus bhe package has been 
wxibben such bhab bhe user does nob have bo call upon special 
pr ogx'ajnmlng skills or knowledge^ excepb fox kriowing elemenbary 
concepbs aboub operabiori of a PC and giving execubion commands 

The package pX'Ogrammes liave been developed on an IBM PC 
compabible The px'Ogx-amme is wribben in PASCAL The graphic 
visuals and Inbex'acbion are accomplished by TURBO PASCAL - vei^slon 
4 0 The edlblng facillbtes for inpub/oubpub is px'ovlded by giving 
access bo WORDSTAR 2000 from wibhln bhe package ab requli^ed places 
and re-enberlng inbo bhe package ab bhe same place The currenb 
vex'siori of bhe package and ibs subsequenb modlf Icabions are 
expecbed bo x^un on IBM PG-XT or compabible wlbh oi' wibhoub a 
colour monibor The package ibself is pr-obecbed fr^om unaubhorlsed 
use by a password code wibhoub wlilch ib can nob be opened by bhe 
ixser 


The algor Ibhm of bhe package pr ogr amme can be besb 
represenbed by a macr-o flow char^t presenbed iri Pig 4 1 The logic 
can be bioadly divided inbo Cl> package access code^ cover* page 
and geriex*ai ublllbies^ <2> menu genex ablon and problem 
idenbif icabion Cselecbion of modelx^&lmulabion bechnique and oubpub 
foxm> from bhem^ C2'> accepbance ol inpub^ display arid edibing of 
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Input. ^ and checking at input. ^ <4> execut.ion of appx’opr4at.o 

px'O^namme^ <[5>dis;pl£iy of dl^ilaL^^r-aphlcai out^put.^ <6> <&ri- queueing 

Ihe outrput.^ and C7> saving of out.put. on a file in t,he user* amea 
fon fut^una nerer-ance and T*at.nieval 

For* speedy pnocessing: arid lo avoid any memor-y p3>oblems^ 
Ihe eriLine package ppo^X'amme Is made up of dirfenent. pr-o^a^ammes 
wWch ax*e execubed wibh pr*opei> cont.x*oLs t.hz*ou^h a babch file in MS 
DOS environment. CMS DOS version 3 3 or Latrer> 

4.2 Menus: and Problem Ider^tiricatlon 

All seleobable items at a menu are written in individual 

windows and the menu ldentlf*loattori Itself occupies a separate 
window The selection of a menu item is effected through an arrow 
^ "^4^^ key on the keyboard The forwar*d movement of the user 

through the list of menus is affected by successive selection of 

the desired menu Item^ the backward movement of the user is caused 
by pressing P2 or through selection of STAY option The package is 
designed to provide the user ’’on- line- he Ip “ to ^ulde him tlirou^h 

the package and a ^ulde option to assist in problem 

identification The help is obtained by selecting HELP or px'essln^ 
Pi while GUIDE Is a menu selectable item The other options 
provided for easy movement in the packag^e include NEXT CF3> option 
to continue with curr ent selection^ QUIT CF4> option to quit or 
exit from the package arid PRINT CF5> option to px*lnt screen 
display 

4 3 Selection Identification and its Implementation 

Each terminal selection Is Identified by a rmlque 
identification number and It^s referred to by that number 
throughout the package px*o^ramme execution The various terminal 
selections of the package and theix* identification numbers are 
^iven in Table 4 1 Each identified problem type is programmed 
according to the concepts and pr*oceduTes described in Chapters 5 
and 6 
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Table 4.1. IdenblflGa'tiori Numbers and Dummy Files Assoclaled 
wlt-h Various Terminal Selections. 


Selec'tion 

IdenLific£iLlon 

Number' 

Corresponding 
Dummy File 

UNIS-DIGITAL 

1100 

ALG.ll 

-GRAPHICAL 

1200 

ALC.12 

MULT YS- REGPGONG 

2100 

ALG.21 

-CONTOUR 

2200 

ALC.22 

PUFFS-DIGITAL 

3000 

ALC.30 

-GRAPHICAL 

GlOO 

ALG.31 

SEGMNT-DIGITAL 

4000 

ALG.40 

-GRAPHICAL 

4100 

ALC.41 

PUPSEG- DIGITAL 

5000 

ALC.50 

-GRAPHICAL 

5100 

ALC.51 

GUIDE-UNIS 

6100 

GUIDE.61 

-MULTYS 

6200 

GU1DE.62 

-PUFFS 

6300 

GUIDE .63 

-SEGMNT 

6400 

GU1DE.64 

-PUFSEG 

6500 

GU1DE.65 

EDITOR 

7000 

OPDEM 

DISPLAY-UNIS 

7100 

DIS.71 

-MULTYS 

7200 

DIS.72 

-PUFFS 

7300 

D1S.73 

-SEGMNT 

7400 

DIS.74 

-PUFSEG 

7500 

DIS.75 


t 



■4.4 SimuLatlort TectxnlqueEs 


In adl five slmula-tion t.echrilques ident-lfied £ia "UNIS'% 
•'MIJLTYS*% "PUFFS"*, ""SRUMOT" AND "PUGSEG" employed UNIS 

simple <3eiussien model Lo simulet,e poUutent dispex»sei Trom single 
si^eedy erid contrlriuous. emlsslori souur'ce ’wf'iichi mey b>e e g^i^oiirid 02 > 
eleve'ted point, sounce in e st^et^ioneny end homo^eneotis et.mosphene 
MUl.TYS uses seme model es UNIS t.he modlt^lGet.ion t.Ket. t.iie 

ef^f*eGt. due t,o mui triple sounoes Is considered by lineer 

superpoeiblon of" t.he ef'fecb due bo severed single iridependentr 
sources Furbhex’ debells ex'e px^esenbed in Ghepbex^ 4 

PUFFS, SEGMNT exid PUFSEG ere besed on exbended Geussien 
epproech end ere used bo slmulebe pollubenb dispersel from single 
urisfbeedy ^i^ound or elevebed poinb source in non~sbeblonery end 
Inhomogeneous ebmosphere in predominenbly celmwlnd, brerispox'b end 
li^hbwind condibioris respecbively Mebhemeblcel debeils end 

el^orlbhims for bhese simulebion becHnlcjues ere presenbed in 
Ghepbex’' d 


4,5 The Peckege Progre3f¥xmes 

The pecke^e comprises of* severe! pX’O^remmes vhdch in 
burn consist of* severel procedures end buncbioris Following is e 
lisb of" pro^remmes end bheix- buncbions in brief* 

1 INTRODUG EXE Irdbiellzes peckege by displeyin^ Tirsb 

cover pe^e end deletes sever*el dummy files <^whlch ex^e used es 
control files end cx eeted during the execution ol the pecke^e>^ if 
exist 


2 PERMIT EXE* Dlspleya passcode option end ellovs access 

to the package only if cor^recb passcode is ^iven Else makes the 
f lie DUMY which leads bo bermixiablon of the package 

3 INFO EXE Gives ^enex'el utilities end information about 
the package 

4 MAINMENU.EXE Displays menu of verdous levels^ allows 

selection of a particular' option through fox'ward or backword 
movement^ ^ives help messages to have desired movement in the 
package end guides in making appropriate selections Exit from 
this prog^ramme creates dummy files HLPDUMY, DUMY end MENU OUT 
depending upon the selections made 

5 RUNN EXE Reads various identification numbers from 

MENU OUT end creates corr-espondin^ dummy files CTable 4 1> 
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6 DOLl EXE D^jle'Les DMLVl sind DMLVll^ iT exist. Checks "the 

existence of the p3>opei> Input file if the option selected is 
UNIS— DIGITAL end the Input is thr oug^h e file CINPlDj^ cr*eetes dummy 

file DMLVl ^ if pnopen input file does not exist and is to be 

cneated/edited^ else ci^eates ariother* dummy file DMLVll leading to 

ignonin^ the cunnent selection 

7 DOL2 EXE" Same as DOLl EXE hut fox» selection 

UNIS— GRAPHICAL Files DMLVl and DMLVll ax*e neplaced by DMLV2 and 
DMLV12 nespectlvely 

8 DOL3 E3fE* Same as DOLl EXE but fox* selection 

MULTYS-REGPCONG Files INPl, DMLVl and DMLVll ax*e nepiaced by 
UVP2^ DMLV3 and DMLV18 i*espectively 

9 D0L4 EXE’ Same as DOLl EXE but for* selection 

PUFFS-DIGITAL Files INPl, DMLVl and DMLVll ar*e replaced by INP3, 

DMLV4 arid DMLV14 respectively 

10 DOL5 EXE Same as DOLl EXE but for* selection SEGMNT- 

DIGITAL Files INPl, DMLVl and DMI Vll are replaced by INP3, DMLV5 

and DMLV15 respectively 

11 D0L6 EXE* Same as DOLl EXE but for selection PUFSEC- 

DIGITAL Files INPl, DMLVl and DMLVll are replaced by INP3, DMLV6 

and DMLV16 respectively 

12 DOLL EXE Vlien the EDITOR option is selected;, through 

this pro^X'amme a dummy file DUMl is created;, exlstance of which 
sends the control to the main menu after editing facilities 

13 EDIT EXE* Displays the EDIT-GUIDE and allows excess to 

WORDSTAR- 2000 for editing 

14 GALTABLE EXE* Allows display of stability classification 

<TABLE 1> and coefficients employed for computing dispersion 

standard deviatloris usln^ power laws CTable 2> 


15 GUIDl EXE* Displays ^ulde lines for selection UNIS 

16 GUID2 EXE Displays ^uide lines for selection MULTYS 

17 GUID3 EXE Displays ^ulde lines for selection PUFFS 

18 GUID4 EXE* Displays ^uide lines for selection SEGMNT 

19 GUID5 EXE* Displays ^uide lines for selection PUFSEG 

20 check: EXE* Deletes MENU-OUT and DUM, it exists Allows 
more selection if required by creatiri^ dummy fll^ DUM 

21 LAST EXE It deletes all dummy flies created during the 
execution of the package and displays end cover pa^es 
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4.6. AoxilliM'y £i]nid Fll^ss 

The execut>lon ot t^he package involves Lhe uise oJT 
auLxlllaiivy riles GLBUNIT, UNITDGTL, UNITORP, UNITMUL, UNITPUP, 
UNITSEQ arid UNTTPFSG alon^ wirh soft^war^e MS DOS Cvension S3 oi^ 
larex^>, TURBO PASCAL <vension 4 0:> and WORDSTAR 2000 Each of Lhe 
pnog^nammes mentioned in section 4 5 uses some on all of Lhe 
auxiliar'y files The auxiliary file QLBUNIT contains a number* of 
pnocedux'es and functions to Cly ^ive vanlous types of audio 
sigrials^ C2> display dlffer-ent messages^ CS> permit accepting real 
value ordy^ <4> permit accepting Interger value only^ <[5> accept 
response YES or NO^ C6y read and display stability classlf icatlon 

arid dispersion stairidar d deviation coefficients tables arid C7 > 
check if the desired files exist in the user area/disk Function 
and description of auxillai^y files UNITDGTL, UNITGRP, and UNITMUL 
are presented in Chapter 5 while the same for UNITPUF^ UNITSEG and 
UNITPFSG are given in Chapter 6 

In addition to auxiliary files there are several input and 

output files which are X’-equlredx'created dur ing the execution of 
the package These files can be classified in two categories^ <!> 
dummy files and C[2> real files Dummy files are created to 

trarisfer proper control from one programme to the other and are 
deleted automatically upon the termination of the package ox» in 

between Real files cari be accessed by the user and can be read 
from outside the package Names of real flies can be supplied by 
the user 
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5 UNIS AND MULTYS 

POLLUTANT DISPERSION IN STATIONARY AND HOMOGENEOUS ATMOSPHERE 

5.1 Scope 

This cKapt^er- pr-esen'ts al^or'l-thms employed in developing 
comput^ei> prognemines: fon t,he slmulA'tlon or concen'tnet^lon 

comput^fii'tion or g£ifi:eous poliurjMirs x-eie&sed r>om ®-teady 
(elngle/'multiple^ gnoundL/elevated^ point, soimceCfi:> in fct^ationany 
end homogeneous atmosphere using the basic Gaussian approach A 
brier description or the concepts and mathematics involved is 
given to assist in understanding the algorithms 

5.2 Gaussian Formulation 

The Gaussian rormulatlon tor atmospheric dispersion 
attempts to simulate the gross behavior or plumes emitted rrom 
ground/elevated sources For localised point sources such as a 
stacks the general appearance or the plxime might be represented by 
the schematic shown in Fig 5 1 Although the plume originates at a 
stack height it rises an additional height AH CPlume rlse>^ 

owing to the buoyancy and momentumn or the gases leaving the stack 
Consequently^ tor practical purposes the plume appears to be 
originating From a point source at an equivalent stack height H 
Corrective stack height) » H^+AH 

2 





Fig. 5.1 Wind oriented coordinate systen indicating 
an elevated point source located at (O^OJ 
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The rcr-mulat^ion Is based on -the appr' 0 >d.maLiori bbat, bhe 
conoenbr ablon dlstrX^ibublon downwind of a points sour-ce in bhe 
atmospher’lc boundax'y layer* Is Qausslari bub wlbh unequal dlspei^sion 
par ame bars in bbe horizonbal <cros&wind> and venbical dlreobions 
For developmenb of" mabhemablcal model a wind orienbed coordlriabe 
sysbem Is employed in wMch bhe x axis coincides wibh bhe 
prevalliri^ wind dlrecblon <!P1^ 5 1^ The ^eomebrlc facbor bh.ab 
represenbs bhe rablo ot bhe concenbrablon ab ari ax'blbrax’y obb-axis 
locabion <x^y>zi> bo bhe concenbr-abion ab bhe cenbreline of" bhe 
plume is bherefore ^Iven by 

r ^ 2 

F<:x,y,z) - €!xp-|- - 1 C5 1> 

2<y^ 2& J 

y z 

Qausslari dis-tr-ibut-lon Is used In bhe abmosphei>lc boundeiry Layer- 

- 2 

even bhou^h bhe burbulence bhere is nob isobroplc <e ^ ^ u^ 5^ 
w'"^> oi' homogeneous Ce ^ ^ w^^ is ^ funcblon of Unequal 

laberal arid verbical dispersion coefflcienbs are used bo 
compensabe^ ab leasb in par^b for bhese deparbures 

The ^eomebr Ic facbor F<;x^y^z> has bhe value of unlby 
alori^ bhe plume axis where y * z ■ O and bherefore represenbs bhe 
qxranblby whdch should be mulbipUed by bhe axial concenbiablon 
c<x^O^O> bo yield bhe concenbr ablon ab any arbibrary poslblon 

c<x^y^z> m cCx^O^O> F<x^y^z> CS 2> 

The axial concenbrabion c<x^O^O> Is debermlned by considering 

species Gonservabion i e source sbrengbh 0 Is equal bo ibs flow 

X abe bhrou^h any plane normal bo bhe x axis This requires bhe 

assximpblon of sbeady and conbinuous emisson The species flow rabe 

through ari elemenbal ar^ea dA normal bo u is cudA and bherefore bhe 

conservabion of mass x-equlres bhab bhe flow bhx-oug^h any downwind x 

plane is equal bo bhe source sbrengbh 

+00 +00 

0 ■« u cCx^O^O / X FCx^y^z^ dy d’X CS 3> 

“00 -00 

The Inbe^i al is readily evaluabed as 

+ 00 +- 0 D 

/ / F<x^y.z> dy dz *■ 2 n o* <xl <y Cx.} 

' j'j y SZ 


- OD ~0O 


C5 4> 
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Thus hhe cen-tx^ellne concerihr*ahlon can be given by 


cCx^O^O 


Q 


2n UL <y (y 

y z 


C5 5> 


The spablai vaidabion of concenbi ablon downwind of bhe sounce cari 
be obtained by Insenblng Eqs 5 1 and 5 5 info Eq 5 2 


cCx^y^z:) 


2n u o o 

y 2 


exp 


{- 




3C 


C5 6> 


Eq 5 6 gives "the basic Gaussian plume formula which can be used 
bo compube downwind concenbi abion fx'om an isoiabed poinb source 
locabed ab bhe origin in a medium oi infinibe exbenb if 

bhe values of dispersion coefficienbs are properly chosen The 
concenbx abion fox^ an elevabed poinb soux^ce can be f ound by 
relocabing bhe soux^ce ab bhe poinb CO^O^H> and by a branslabion of 
bhe Goordinabe sysbem where ^ is x-eplaced by as-H 


cfx^y^z!) ■> 


2n u <y cr 
y 2 


exp 


{- 


20* 




2<y 


-} 


C5 7> 


For a ground level source^ bhe verbical spread of bhe plxime is 
resbricbed in 0 bo region arid hence bhe downwind concenbrabion 

expression bakes bhe form 


cCx^y 


Q 

rr u o O' 

y z 


exp 


{- 
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CS 8> 


Eq 5 7 for an elevabed source is of ben expanded wibh a r^eflecbion 
berm if bhe sux'face is nob a 100% sink for poilubanb The downwind 
concenbr* abion wibh reflecbion is debei'mined mabhemabicaliy by 
linear' super posibion of bwo Gaussian bype concenbrabion curves^ 
one cenbxed ab H Creal> and bhe obher ab -H <limaginar*y> The 
sbrengbh of imaginax'y soux'ce is baken as some fracbion ^ 

parameber xelabed bo bhe chax'acberisbics of bhe ground surface and 
poUubaxib^ of the sbrengbh of real sour-ce Eq 5 7 with reflection 
berm bakes bhe form 





2fT u o* 

y z 






2 

y 

Cz-H>^ 

2 

2 1 

20* 

20* ^ 

y 

z 

2 

2 


Cz4-H> n 

1 2 

2 I 

^ 2<y 

2iiy J 


5.3 Appllcjaitlori at G^xiscssi&rs. For^molA^ 


cs c» 


Applicfa'Llon at -Lhe fiif‘cr*ement.ioned GaiussiaLn 
x«equlr«es inToi'nriaition about, pollubarib eml&sion r»at.e 0^ offeobivo 
sback height, wind veloolby disfporslon paramebar-s and <y^p 

and neflecblon p£U'amat.er P ^ While emisrslon nat^e can be 

l>€it 

esbimaled/measuned iTalrdy accux-abely^ esblmatlons/meaisux’emenbs: of* 
<y^ and have bo be based on some assumpbions and 

Ideailzablons The following aecbloris descr-lbe bhe 

techrilquey'pnocedune adopbed for* reasonable esbimabes of bhese 
pax^ameber s 


Effecblve sback hel^hb Effecblve hei^hb H of bhe vli>bual 
oni^in of plume i& obbained by adding a benm AH^ bhe plume nlse^ 
bo bhe acbual sback hel^hb Thene ane numer-ous mebhods fox* 

caiculabln^ AH^ and bhese ane discussed ab some len^bh by Sber*n 
<;iC>68> Basically^ bhr-ee sebs of pax*ameber*s conbnol bhe nise of a 
gaseous plume Injecbed Inbo bhe abmosphene fnom a sback These ar-e 
sback cha3>aGber*isblcs^ mebeonolo^lcal condibioris^^ and bhe physical 
and chemical nabur-e of bhe effluenb A lar*^e number* of analybical 
expnessloris have been proposed bo r*elabe bhese facbons bo plume 
X’lse predlcblons Carson and Moses compared 711 observed 

plume rise values wlbh calculabed values ^Iven by eleven dlfferenb 
equabloris Their* conclusion was bhab bhe following bhree equablons 
fox* unsbable^ neubral and sbable abmosphere predlcb reasonably 
accurabe AH values 


AH 


AH 


AH 


V D S 

C ri J 

3 47 4- 5 15 <unsbable> 


u 


u 


V D 

o 95 ^ + 2 64 — Cneutr>ail> 


u 


u 


V D Q O 5 

-1 04 ^ — 4- 2 24 ^ <]sbable^ 


u 


u 


<5 10> 

<5 11> 

<5 12> 


Here is bhe sback ^as exlb veloclby in m^s^ Is bhe sback 



Zi 

exit dlfuneler in u in m/s £irid Q. is 'the the heat emission r»ate 

s h 

in Kjys 

Wind velocity The wind velocity varies with altitude Czy and 
the typical variation or u with !z can be obtained by following 
equation which includes a factor to adjust for various stability 
conditions in the atmosphere 



C5 13> 


Here u and u are the wind velocities at heights z and 
respectively^ and R is the stability pai^ameter which is i elated to 
stability condition of tht? atmosphere CWar k %iid Warner ^ ly81> The 

appropriat.e value of u to use in the dispersiori t:?quatlons Is the 
mean value taken through the plume CTurner ^ 1969 > In most cases 

it would be impossible to determine the mean^ since sufficient 
atmospheric data would not be available In lieu of thls^ the 
average wind velocity at the top of the stack is commonly \ased 
Since in most cases not even this value is known^ the measured 
meteorological value at 10m is used in conjurictlon with Eq 5 13 
to estimate the wind velocity at the plume height 


Dispersion standard devlatloris* As might be anticipated from 
the physical description of the diffusion problem^ the horizontal 
and vertical dlspex'slon standard devlatioris^ <y^ and are a 

furictlon of the donwlnd position x as well as the atmospheric 
stability conditions Many experimental measurements In the 
atmosphere have led to an evaluation and cor relation of <y^ and 
values There axe several sets of charts for these two parameters^ 
and the range of stability condltloris covered in the different 
sets do not rior mally coincide 


Dispersion 

standar d 

deviations are 

determined by 

the 

mixing or dispersive 

power 

at 

t,he 

turbulent flow within 

the 

atmospheric boundary 

layer* For 

t-hls 

pui-pose 

bhe st-a'te of 

the 

boundary layer is cliar’acterized 

by 

P£is quill 

ClQ74:>f as beln^ 

in 

one of the nine 

A-B, B, 

B-G 


G-D, 

F> classes which 
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r^epoi-t^d in text books Ce*^- Wark and Wax^ner^ 1081^ 

Dobbins^ 107^!>. Determinants of' the dlspex-sion classes ax^e wind 
speed arid the state ot incoming ox* out^oin^ solar radiation which 
attect the lapse rate^ presence ox* absence of convective activity^ 
and the dyriamlcs of" the nrdxed layer. 

Fox’ computational pux'poses it is desirable to have the 
two standax’d deviations expressed in mathematical for^m. Due to 
logarithmic riatux*e of" the <y variations with downwind distarice x^ a 
reasonable cxix’ve fit is obtained by power* law expressions. Thus 
algebraic representation of and is frequently expressed in 

terms of the expressions of the type 


P c 

o* » a X and o* ® b x^ dx. <y and O' in mi 
y !z ^ y a ^ 

196P> 


Tadmox* arid Clur^ 


or 


a X and o 


c X f <lx in km> o and o in m; 

y 'z 

Martin, 1P76> <:5.15> 


where a, b, c, d, f, p and q are empirical constants which depend 
on the stability class and wind velocity <!Wark and Warner, 1^81 1 
Dobbins, 107P>. 

Reflection pax*ameter**. As mentioned eax*ller, reflection 

parameter* P . is Influencd by the cliaractexi sties of the surface 
rei 

arid the pollutant. No general guidelines ax*e available for 

selecting the value of P ^ and has to be based on intuition and 

ref 

experience. 


5.4 Option UNIS 

This option is for predicting the downwind, cross wind 
and vertical cocentratlon profiles due to steady arid continuous 
emissions fx’om a single point source under* stationary arid 
homogeneous atmospheric conditions. It has two sub options DIGITAL 
and GRAPHICAL to give output in digital and graphical form 
respectively. Correspondingly thex*e are two separate programmes 
UNISDGTL and UNISGRP. UNISDGTL uses units GLBUNIT <Ref. Sec.4.6> 
and UNITDGTL CFlg. 5.2> while UNISGRF uses GLBUNIT and UNITGRF 
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<F1^5 3> Urdtis GLBUNIT, UNITDGTL aind UNITGRF consist of several 
sLibpro^X’ajYones Fl^ 5 4 presents the schematics of the lo^lc 
adopted In UNISDGTL and UNISGRF which constitutes option “‘UNITS" 
Details are availahle elsewhere CNema and Taxe Nema^ 19B7> 

FoUowin^ Is the bxdef descrdptlonx^f unction of vaxious 
subprogrammes used in units UNITDGTL and UNITGRF 

Cl> INPU1\ INPOOT* All data are read through these 

subprogrammes and can be ^Iven through an input tile or terminal 
in an inter active mode arid can be stored in a file depending upon 
usei’S choice The input data consists ot the following 

a Source - emission X'ate^ heat content of the gaseous 

stream^ e>d.t velocity^ physical stack height and stack diameter 

b Meteorology - wind velocity measux^ed at a specified 
altitude^ erivlx'onmental lapse xate^ data about day <!:> oi night 

solar- radlatlori 1 e str-ong model ate <.2> ox slight C3> 

and cloud covex^age 1 e clear- Cl> oi- overcast sky c2^ 

G Domain - specifies minimum and maximum distances 

Calongwlth the interval:) in which the pollutaxit concentration 
profiles ar-e deslr-ed 

2 W3^ W2^ W1 These subprogramme s make diffex-ent windows on 
the screen for- desired patterri of screen display 

3 OPENFILEl, OPENFILE2- OPENFILEl assigns and resets the 
input file to INPl while OPENFILE2 checks whether the user wants 
to store the output In a file or not It cari store the output in a 
user specified file 

4 VELOGT* This subprogramme computes wind velocity at any 
altitude accox-dlng to a powei law by selectlxig appx-oprlate 
stability pax-ametex r-elated to the stability condition of the 
atmospher-e 

5 STBGLS This subprogramme selects appx-opilate stability 
index Ci to 6 fox- stability classes A to F respectively > 
cor i-espondlng to the meteox-ologlcal data and wind velocity at 10m 
height 



UNISGRF 
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Fig. 5.3 Progranne 
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ISTOP 


Fig, 5.4 Schematic representation of the logic adopte 
In developing option ^UNIS^ 
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6 OLASS^ GLASSl* Subpi’C^r-amme CLASS decod^jfis Lhe S'Lfibili'ty 
class and Lhe subpno^namme GLASSl wxdLcs the same in a ouLpuL 
file 


7 SlGf, SICMA The vardous empinlcal consLanLs involved in 
power- law fox- compulln^ dlsper-slon sLandand dviafioris ax-e assigned 
in SIGMA wldle in SIG^ and ax-e compuLed 

8 BOUYNC* This subpx-ogx-amnie compuLes Lhe plume x-lse arid 
effective sLack height, Appx-opx-latre expr-ession for- plume 3 >lse is 
selecLed depending upon the stability class Stability classes A 
to G are considered as uristable^ D is consider-ed as neutral^ and E 
and F are considered as stable 

9 NOTREF, REFLN, GONG Subprogrammes NOTREF and REPLN 

compute the crosswind and ver-tical conceutr-ation profiles in the 
specified domain at eqpxal Interval for a pollutant released from 
an elevated point sour-ce Goncentr-ation is computed with REFLN 
when the plume implriges on the surface The n^axlmum crosswind 

distance is selected corresponding to 10% of the plume centre line 
concentration on either- side subject to its not being gx-eater 

tliari 2 0 km Signlflcarit ver-tlcal distarice is taken as the 

distance between plume centreline and the ground levels on either 
side of the plume centr-ellne at a particulai* downwind distance 

The subprogramme CONGO is similar to NOTREF except Cl> 
wind velocity is taken at Sm height^ C2!> vex-tlcal concentx-atlon 

pr-ofile is computed fr-om surface level to a distance where 10% of 
the sux-face level concentration is obtained subject to its not 
being greater than 2 O km To avoid floating underflow ex-r-ors arid 
save on computational effox-ts^ the concentration is set to zero if 
the exponent is less than -25 in all the subprogx-ammes. 

FORMl^ FORM2^ FORMll^ FORM12 These subprogrammes are 
included to write crosswirid and ver-tlcai concentx-atlon profiles on 
terminal and output file in digital form 

FORMGl^ FORMG2 These subprogrammes display on the screen 
crosswind and vex-tical coxicentratlon pr-oflles in gx-aphlcal form 
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5 5 Optiom MULTYS 

TMs optxlon Is for predicting pollutant concentration in 
a stationax'y and homogeneous atmosphere usin^ Gaussian approach 
CSectioris 5 2 and 5 33 for a pollutant continuously released from 
multiple point sources C^round level or elevated!) at constant 
emission rates It has two sub options REGPGONO and CONTOUR to 
^ive pollutant concentratioris at selected locatioris in digital 
form and isoconcentration lines in the selected donrain 
respectively However^ the second option is not completely 
developed The first option is Implemented tlirough a programme 
named RECPGONG and uses units GLBUNIT and UNTTMUL which in turn 
corisist of several subpx’Ogr ammes <Plg 5 53 These subprogrammes 
are similar to those used in UNISDQTL/UNISGRF except for some 
minor chaxiges as follows 

Cl> INPUT^ INPOOT “ Source characteristics are to be given 
aiongwlth location for more than one sources^ and receptor number 
and locatioris ar^e to be specified instead of domain 

C23 CONCH ~ This computes concentration due to elevated 
source considering both reflection and no reflection 

C33 SIGMA ~ Subprogrammes SIG and SIGMA <Fig 5 23 are 
combined Into this subprogramme 

<C4> Goncentration contribution from a source is computed only 
If the receptor is located downwind of a particular source The 
aggregate concentration at a receptor is computed by linear^ 
super position^ neglecting the interference due to multiple 
sources 


Schematics of the logic adopted In developing the 
progr amme REGPCONC is presented in Fig 5 6 



RECPCONC 
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iSTXRT 

i 

INPUT : METEOROLOGICAL PARAMETERS 
NUMBER OF SOURCES AND THEIR CHARACTERISTICS 
STABIUTY CLASSIFICATION (TABLE 1) AND 
DISPERSION STANDARD DEVIATION COEFFICIENTS (TABLE 2) 
DOWNWIND TERRAIN CHARACTERISTTCS 
NUMBER OF RECEPTORS AND THEIR LOCATIONS 

, I 

t^RINT TABLE 1 AND TABLE 2 IF REQUIRED! 

I 

COMPUTE WIND VELOCITY AT 10m AND 
IDENTFY STABIUTY CLASS 

1 


SELECT A RECEPTOR , 



STOP! 


Fig 5 6 Schenaiilc represen'ta'tlon of fhe logic aolopfed 

in programme 'RECPCDNC' 
















6 PUFFS, SEGMNT AND PUFSEG POLLUTANT DISPERSION IN 
NON-STATIONARY AND INHOMOGENEOUS ATMOSPHERE 


6.1 SIcope 

The ul'tlmeihe of hhe dlspex slon modelling is "to 

px’ovlde a met^hodolo^y ho predict poUutairit diapersal in a real 
life situation This involves complete axialysis of dispex sion 
phenomena in non-statlonax y and inhomogeneous atmosphere for a 
pollutant leleased from a non-steady source dex’-tainiy^ a problem 
of this natux^e caxinot be analy'zed without simplified assumptions 
and Ideallzatloris In this study atmosphexdc dispex slon phenomena 
for aforementioned conditions is simulated usin^ e>itended Gaussian 
appr»oach This chaptex' briefly describes the concepts and 
nxathematics Involved in developing computer' pX'O^X'ams t or the three 
f ox'mulations namely PUFF <lox essentially calmwlnd condltlons>^ 
SEGMENT Cfor transport conditions^ and mixed PUFF-SEGMENT Cfor 
llg^htwlnd conditions^ 

6 2 PUFF FormuLatloin 

The basic Gausslari approach CEqs 5 7 and 5 described 
in previous chiapter ^ives a plume formula whose validity requires 
the main assumptions of <!:> spatial homogeneity^ Cil> stationary 

conditions and Clil> steady continuous emissions Fux'ther^ these 
equatloris cannot be used dlx'ectly in calmwind condltlorus where j u j 
m 0 However^ simplicity of the Gaussian appioadv its X'elatively 
easy use with clearly measureable meteoi'olo^ical pax'ameters and 
especially the elevation of this method to the quantitative 

decislon-Gonti oiling level have stimulated i esearch aimed at 

removing the limitations of the Gaussian theory to treat the 

complex situations of the real world In pax'ticuiax'^ the PUFF 
approach CLamb arid Nelbux'^ex ^ 197l:> has been successfully applied 

to pseudo steady-state conditions Several studies cZannetti^ 

Shelh^ 1978; Ludwig ^ ert al ^ 1977> have discussed in detail 
the PUFF modelling appx'oaclu impx'ovln^ its application features 
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Th€^ most, obvious way oi" t,r'eat,lri^ non~st,atloiri£o>y 
Gondibioris in emissionsymeleoi-olo^y^ Is bo descr^ibe bhe point, 
sour-ce emission inbo bhe abmosphene by a sei>iea of *‘irisbaribaneous*’ 
puf'f’s whose ciianacbenisbios ane updabed ab each dlspension 
Inber-val Ab Mebeox'olo^icai bhi^ee dimensional fields <wind and 
bui»bulenb sbabus> and emission paxamebei’S ax^e allowed bo change ab 
each mebeonolo^ical bime sbep ^^n\eb^ l)Vi"i£imics of" each puff 
Gonslsb of Cly ^enex-abion ab bhe sour-ce^ C2y plume nise^ <3> 
bnarxsponb by advecbive wind; <i4> diffusion by abmospher-io 

bui’bulerioe^ C5y ^nound deposibion^ dny and web^ and <6^ chemical 

bnansfonraabion Each puff is chai^acbei^ized by bhe following bime 
vax^ying pai ameber-s evaluabed ab ibs cenbnal poinb 

C^CPX^ GPY^ GPZ> coox'dinabes of bhe cenbre poinb of puff 


h 

P 

M 


elevabion of cenbne poinb above gi>ound 

tin flab benxain h » GPZ> 

P 

pollubanb mass in Puff 


o, 


Y" 


sbandard deviabions of bhe Gaussian concen- 
bx'abion disbribubion 


Generablori of plume^ Ab each bime inberval Ab^ a new puff is 
added bo bhe puff ^chaln^ fnom bhe sounce The par^amebexa defining 
each new puff tiave bhe following iniblal values The cenbx'al poinb 
Is seb ab bhe souxce^s exlb poinb plus bhe vexd leal plume rise AH> 
M w 0<b> Ab^ where 0<^b> is bhe eux^x^enb emission x^abe^ and o ^ 
I'epTesenb irdblal o^s of bhe plume ^ ^nd <y^ « 

AH/3 16, Zannebbi, 1986 > 

Adveebion Ab each bime inberval Ab^ cenbi^al poinb of each 
exlsblng puff is adveebed according bo bhe Gux*x«enb wind veebor 6 w 
tu,v,w> ab puff^s locablon as follows 

CGPX,CPY,GPZ>^'^^'^^ - CGPX.GPY^GPZ)^*^^^^ + 3 At- C6 1> 


Diffusion During each Ab bhe puff^s o^s are basically 
Gompubed usixig power law ejcpressioris tEqs 5 14 or 5 15> in which 
puff^s dix’ecb dlsbarice fx'om bhe soux’ce is subsbibubed for x 
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However^ > vlr-t^ual distance <t,r»arispor*i, condit,ioris:>/vla:'t,u£il 
C^calmwlnd coridlt,loris> concept, is used whenevei* met^eoi olog^y <i e 
stability class:) is differenl at pufl'^s new and old centne point 
either due to non- stationary or non- homogeneous condition of" the 
atmosphere CLudwi^ ef al ^ 1977 ^ Zannetti^ 19Q11 


Virtual distance concept for <y computation* 1 Current o 
fxiriction <corr espondin^ to new stability ciass3 is equated to 
^ ([old> calculate appropriate virtual dlstarice Cd for o and 


H 


d,, for o ) 
V !Z 

!Z 


Cold> 1/q 




r 1 


1 

r 

y 1 

- and 

d ■■ J 

1 2 



1 b 1 

1 J 


z I 

1 J 


in S>I(3MA1 


C6 2> 


Cold>^l/b 


and d 


H 


<:old> ^^l/d 
O' -f ^ 

•Z 


V 


in SIGMA2 


C6 


2 New values of o^s are obtained by substituting <dv + U At> for 
X in Eqs S 14 or 5 15 


Virtual a^c concept for o* computation 


In calmwlnd 


conditions (!i e U < U , . U is taken as U , > o^s are considered 

min min 

as functloris of time <mor e exactly, the a^e> in power law 
expressloris <;Eq S 14 or 5 15> instead ot downwind distance Simlar 

to virtual distances, virtual a^ea Ct for o and t for <y > ar^e 

H ^ z ^ 

computed as follows 

Cold)_l/p ^ 


H 


<y 


<y 


Cold>^l/q 


and 


a U 




in SIGMAl C6 4> 


min 


b U 


min 


a 


<o Idl^l/b 


U 


min 


and 


o* 


<o 1 d :> 


o U 


min 




IZd 


in SIGMA2 


C6S} 



New Values of 




ai»e 


ob'Lairi.ed by l ollowin^ expx'eaaiori 


)r r espoiidiri^ bo At) 

<y Ct.-hAi.> m a + At,:>^ and cr Ct,+A-t> » b Ct,„ + Ab>^ 

y min ss 


Cln SIQMAl:> 


C6 6 > 


or 


o-y<t+At:) - a and o-^Ct.+Ab> - b Cb„ + + 1 

H 


min V 

2 

Cln SIQMA2> 


C 6 7 > 


Dr-y 

and wet 

deposition 

Both dry and wet 

deposition 

for 

the 

pollutant 

ar e s-lmulated 

by 

f irst“-ox der reaction 

schemes 

and 

are 

computed 

during 

each 

At 

by ari exponential 

X'eductlon 

of 

the 


pollut/anb mass 


. .Criew> .^Cold> 

M « M e>^p C~ Pj AbySdO^OOOl C 6 S> 

Here J indlcabes dr*y CJ*«1:> oi' webCJ*2> deposlbion and is bhe 

cor r espondin^ percentage of reduction per hour* C%y'Kn> P and P 

X ^ 

<^ari be obtained fnom deposition velocity values and px ecipltation 
datii as follows 


P^ « 360,000 V^x^AZ and ^ P^/CIO T^y C 6 

whene is the cunnent pollutant deposition velocity Cm/s^ at 

element's locatlori, A 2 is the vertical thickness of the puff CAZ 
Bi 2 ^ the pollutant scaven^ln^ ratio, P^ Is the cur-nent 

pi»eclpltatlon nate at element's location Cmm/laO and Is the 

current thickness Cm> of the precipitation layer at element's 
location Dr^y deposition is consldex'ed only when the plume has 
reached the ground Cl e 2 o > H“ZR:> 

Ghejridcal transf ox-mation Durdn^ each At, a flx^^st- cruder 
chemical reaction scheme is adopted, in which the chemical 
tr ansf ox’matlon term r*educes the mass of pollutant in each puff 
accox'din^ to 


j^Cnew> ^ exp C- K Ab/3A0,000> <6 10> 

where K is the poILutarit chemical tr ansf ormatlon factor' expressed 
as a pex'centa^e of reduction per hour' C%/hr^ 
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The AG value is compu'ted ordy if it is significant Ci e > 0 001 

This is ensux-ed by compaiing the i-atlo ^computed 

3 , 

using Eq 6 11 by subetituting AG-0 001 pg/m , <y^ ■ o-^f and 

appropriate puff and x-eceptoi cooi'dlnates> and CD/'o')^ <lx'atlo of 
t-he direct distance between puff and receptor to or o■_^^ 

whichever Is greater > If CDy<y>^ > AC is set to zero Eq 

All is expanded with a reflection term Cslmllar to that in Eq 

whenever 2 0 - > CP2-ZR 

z 

Puff splitting Breaking of a plume into pufls allows the 

evaluation of their dynamics as a furiction ol local time-varying 
meteorological conditiorus In particular > during each At, the 

centre point of each puff moves from an old to a new position 

<advective dlsplacement> Large advectlve displacement due to an 
Increase in wind speed may affect puff's abilty to represent the 

continuous plume by reducing resolution The splitting technique 
<Zarinetti, 1981), which generates, when required, a sufficient 
number of fictitious puffs along the puff's trajectory during At 
to maintain sufficient resolution, is used This splitting of puff 

Is performed for computing AG when Cl> AG>0, and C2> advectlve 
displacement >15 0 - In this splitting computation mass M of the 
pulf is equally distributed among the split puffs along the 

trajectory from the old position to the new one and values 

for the fictitious puffs are irxterpolated linearly between the 
respective o values at the puff's old and new centre point 

Computer programme using PUFF formulation Main objective of 
this progr-amme Is to predict time dependent point concentration of 
a pliutar.t released from non- steady source under essentially 



Jt) 

calmwind conditions The pnognamme named "PUPFS" consist of 
sevenal subpr-ognammes CFig 6 !:> and its lo^-ic Is schematically 
presented in Fi^ 6 2 Following are some salient featxxres of this 

pr'0^r’£tmm€* 

1 Erit/li'^ r-eg^iciri of soul}>c<& Irif’luerice is dividi&d lnt,o 
mei^eoi oio^ic^d ^x*ld corislderdn^ the sour^ce to be situated at the 
ceriti»e of the left ed^e of the middle X’ow Henoe^ ther'e will 
always be odd number^ of 3>ows In the ^nld Each block of the ^nid 
is allowed to take dlffer-ent meteox^olo^ical parameter’s after* every 
time step 

2 The plume profile specifies puff number^ its centre point 
coordinates > pollutarit maiss contained in it and the corresponding- 
dispersion standard devlatioris 

3 Loss of pollutant duo to chemical reaction oi dx’y/wet 
deposition Is computed assuming a fix’st ox dei x-eactlon scheme at 
4-voi’y pulf genex’atlon time step 

4 Concentration computations are carried out at user 
specified time Interval This implies that the siibprogramme CONGTR 
is called only when concentration values are I’equired 

5 To avoid excessive computational efforts^ maixlmum number 
of fictitious puffs Is X’estricted to 50 Further to check floating 
underflow error’s^ exponential tex’ms are set to zero whenever the 
exponent becomes less than -25 

Subprogrammes VELOGT, STABGLS, BOUYNC, W3, W2, W1 and 
OPENFILE2 axe same as ixsed In UNISDCTL OR UNISGRF A brief 
descr Iptlonyf unction of other subprogrammes is as follows 

Cl> INPUT The data read through input Hie <:iNP3> ax’e as 

f oilows 

a Source - same as in "UNISDGTL** 

b Meteorology — In addition to the data read in 
*‘UNISDGTL*% other data those are to be read are the angle of wind 
velocity vector in tlie vertical and horizontal plane from x axis^ 
deposition index Cl for dry and 2 for wet>^ the rate of 
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1 

INPUT : METEOROLOGrCAL I^AMETERS, SOURCE 
CHARACTERISTICS,STABILITY CLASSIFIGATK>N(TABLE1) AND 
DISPERSION STANDARD DEVIATION COEFFIGENrS(TABLE2) 
DOWNWIND TERRAIN CHARACTERISTICS & EXTENTS, NUMBER 

OF RECEPTORS AND THEIR LOCATIONS, POLLUTANT 
CHARACTERISTICS AND STUDV DURATION 
I 

SELECT A SET OF SPECIFED METEOROLOGICAL PARAMETERS! 
FOR INDIVIDUAL BLOCKS OVER ENTIRE GRID 
COVERMG TERRAIN EXTENTS 

I , 

COMPUTE WIND VELOCITY AT 10m HEIGHT AND 
PLUME BOUY^NCY AND EFFECTIVE Hs 

, I 

COMPUTE WIND VELOCITY AT K(5in F Hrf5m) 

PLUME BOUVANCY AND EFFECTIVE Hs 


i 



ADVECT ALL EXISTING PUFFS WITH VELOCITY AT 
CURRENT LOCATION AND GENERATE 

NEW PUFF AT SPECIFIED INTERVAL 




COMPUTE CURRENT PUFF MASS AND CROSSWIND AND 
VERTICAL DISPERSION STANDARD DEVIATIONS 


ARE RECEPTOR CONCENTRATIONS REQUIRED ? r 


r? 

SELECT A receptor! 


.1 

COMPUTE SUM OF CONCENTRATION! 
CONTRIBUTION DUE TO ALL PUFFS 1 

1 

SELEC 
NEXT REG 

-^^RE ALL RECEPTORS CONSIDERED ?1 

EPTOR ~j 


ARE METEOROLOGICAL CONDITIONS CHANGED ? 

yI 

^IS STUDY DURATION OVER ? 

. yJ , 

OUTPUT: PRINT/PLOT OF PLUME PRORLE 
RECEPTOR CONCENTRATION IN TABULAR FORM 

stop] 


Fig. 6,'d Schenatic represen'tatlon of fhe logic adopted 
in programme TUFFS'' 
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preclpltat^ion, pollutant depoeltioii velocity arid thlckne&& of 
pnecipltation layer* Theae data ai’e to be provided foi* every block 
of the gr*ld and for every meteorological time step Instrument 

minimum vind velocity is used as a criteria to determine calmvind 

ooritdi'tlorus 

c Dorrtaln of t,he ^vid - M£od.murn values in t,he x and y 

dir'ect.ioris wi"th *the com espondln^ block sltze 

d Recepbox’S ~ numbei* of" r^ecepbox^s and t,h.eix» coox'dinatrea 

e Time vax^lables — maximum bime at inCluence^ 
meteox'olo^ical bime sbep^ concenbnablon calculabion and puTf 
^enenabion inbenval 

f Obhen - pollubanb scaven^lri^ rabio^ chemical 
br arisf ox'mablon facbox* Tor* bhe pollubanb^ coefflcienb of X'-ef'lecbiori 
f X’om ^x-ouxid 

STBLTY* Ib finds oub bhe sbablliby class fon evex^y block 
of bhe ^3>ld ab evei^y mebeoi>olo^ical bime sbep on bhe baLsis of* 

Pasqulll sbabillby classif icabion 

<3> SIGMA This subpr-o^x^amme is used bo compube dlspex'slon 
sbarulax’d devlabloxis The only modlflcablon ovex* bhab used eanller* 

(ie ^ in RECF'GONO is bhab bhe vlx*bual dlsbaxice <;bx»anspor»b 
coridlblons> vlx^bual fcalmwlnd condlblons> coxicepb is used 

whenever^ mebeonolo^lcal Ci e sbabillby ciass> is diffenexib ab 
puff^s new and old ceribx*e poixib elbher* due bo non-sbabionany or- 
inhomo^eneous condiblons of bhe abmosphex-e fSecbion 6 2> 

C4> ADVEGT Advecblon of each puff duidng a period is done 

bhx'ou^h bills subpno^T amme Wind veloclby vecbon ab bhe px>esenb 

iooabion of bhe puff cenbx-e is compubed and bhen lb is advecbed bo 
bhe new locabion^ hence ib finds oub bhe new coondlnabes and bhe 

block index in which lb Is locabed 

<5> LAST* Ib assi^ris bhe chax^acber-isbics bo bhe newly bonn 
puff, e g bhe coondinabes, dlsper^sion sbaxidand devlabions and 


mass 
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<6> CONCTR’ It- comput-fts t-he aonc^nt-rait-lcin cont.r»ibut.ion t-o 
recept-or* by ^v^ry puTf Only t-hose pufTe cont-nibut-in, 
eigniricaLnt-ly t-o & recept-ox' ane considered If the eff'ectivi 
legrith of the puff is greater* than 1 5 then puff spUttlrig li 

adopted to maintain sufficient resolution, In which the puff li 
assumed to be ar» aggregation of a nmnber of virtual pixffs, theli 
characteristics are computed by linear Interpolation, and pufi 
mass is divided equally among the spUtted puffs Aggregati 
contribution of all the puffs is reported as the flna 
coricentr atlon at a receptor 


6.3 SEGMENT Formulation 

The segment apporoach was basically developed to remove 
limitations of steady-state formulation while maintaining 
simplicity of the Gaussian theory The segmented plume mode 
CBenkley and Bass, lOSO, Chen &t al , i97Q, Hales et <tl , 1^77 
which is dynamic arid considers non-steady state condltioris 
however, still requires transport conditions CU > lmx's> 


Segment method breaks the plume into Independent 
segments whose initial featxrres arwi dyr^imlcs are a function ol 
local time-varying emissions and meteorological conditions 
Therefore, it is able to stimulate non-statlonary anc 
Inhomogeneous dispersion conditions Segments are considered as 
sections of a Gaussian plume Each segment generates b 
concentration field which Is still basically computed by Eqs 5 7 
- 5 9, and represents the contribution of the entire virtual plume 

passing tlirough that segment as Illustrated In Fig 6 3 



S = REAL SOURCE 
s' = VIRTUAL SOURCE 
R = RECEPTOR 

R" = representative OF R .-"-I 



X 
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Thereat only nne clos€?stx> a.ff'ectrS; "tHe 

ccnc€;ri'Li’£i'Llon compii'L^t.lon fat, 4Ef£ich I’^jcepLor* Dyri^amics of" 6?Aoh 

segment, is: siimllfm t^o ttifat, of" each puTf except. t,hat. t.he segment, is 
Ghai>ac’ter'i!zed by bwo points Cleading and tx'ailing> on the 

segmented plume as follows 

1 Goor-dlnates of the leading point GPXaFS:), CPYCIFS^, 

GPZCIFS> 

2 Goondinates of the ti^alllng point* GPXaSSJ), CPYCISS>, 

CPZdSSJ 

3 Standax’d deviations of the Gaussian concenti’ation 

distribution o aFS>, o dSS>, cy aFS> and o- CISS> 
y y ^ 2 2 

The pnocedune adopted for* Ci> the genenatlori of 

segmented plume^ C2y advectlon^ C3> diffusion^ <14> dx^y and wet 

deposition^ and CS> chemical tx-arisf cremation is same as described 
in Sec 6 2 However ^ the cokicentnatiori contidhution due to a 
segment is computed as follows 

Goncentnation computation* The concentration computation due 
to a segmented plume requires identif ication of the segment 
closest to the receptor R arid utilization of the segment's 

variables for^ computings using basically Eqs 5 7 arid the 

concentx'atlon field generated by the equivalent plume passing 

through that segment (!Fig d 3^ The parameter's in Eqs 5.7 and 

5 0 are evaluated in the following way 

1 Segment's var-lables elevatiorO are evaluated at 

the point R" Chase of the prependlcular drawn from the receptor on 
the centre line of the closest segment^ Fig 6 3> Evaluation of 
o and O' is done by linear interpolation Cox extrapolation if R' 
lies outside the first and second point of the closest segments 
using respective 0 'CIFS> and o'CISS> values 

2 O is evaluated as vix'tual cur-rent emission rate i e Q m 

MZAt 

3 Wind velocity is evaluated at R" according to Eq 5 13 

taking atmospheric conditioris at R^ location 



Identrif icA'tiori of Glo&est^ Se^men't The s.e^men'L closjeat, 1,0 'the 
recepl^or* is; Ideri-tifled iri Tollowiri^ steeps 

Sl.ep 1 Distance and of t,he fir*st, aiid second point, of 

the leading segment ar-e computed as follows 

-C<GPXClPS>-XR>^+<GPY<;iPS>-YR>^-KGP2Cir&>- C6 12> 

2 2 2-1 /2 
miicpxass ' y - xk.y +<:cpycis:s2-yr> -kcpz<:iss>-zr> ] Cdi3> 

step 2 If the cux'X erit segmeiit is taken as closest 

s€f^ment ox else is equated to aiid is necomputed using Eq 

6 12 fon the following segment 

Step 3 Step 2 is repeated until the last segment and if no 
s#:^gmeiit satisfies the condition of Step Z, the last segment is 
takeri as the closest segment 

Fvaluation of segment parameters Evaluation ot segment 
j:>a^ ^dimeters for concentr atloiri computations is done fox closest 
segmerit at CEig 6 3> This Is done In following way 

1 Evaluate the effective length CEFL:> and dlnection cosines 
CSl ^ SM and SN) of the centre line of the closest segment 


EFL- r<'GPX<:iF’S>-GPXCISS:»'^-M:GPY<IPS>-GPYCIi>t.>>^ HI- 

<GP2<IFS>-GPZaSS»‘^3^'^^ C6 14> 

SL- <GPXC1FS> - GPX<ISS»/EFL C6 15) 

SM- <GPYCIFS) - GPYaSS)>/EFL C6 16) 

SN- <GPZCIFS) - GP2<IS&)>/EFL C6 17) 


2 EvaluAt-e 'Lhe proJecLlon of t-ho line jolriing second points 
of t-he closes-L segment- t-o t-he i eceptroi- <T> along segment, ctnt-i*e 
line 

T - SL<XR-GPX<:iSS)>+SM<.YR-GPYCISS)>+SN<ZR-GP2<:iSS)> C6 IS) 

3 The Goor-dlri£it,e& of R' ane comput-ed as 

SXR - CPXCISS) + SL* T C6 19) 


SYR - GPYCISS) + SM* T 


C6 20) 





SZR « Gpzass> + SN* T 

4 Th€? <y £ind <y V£ilu€5S. ait^ point aix-e ^Iven by 
y z 


O' a-t R" 

y 


o at R^ 
2 : 


T 

S MO* <CISS> + <0" <!!IFS^ - <y <ISS» =^- 

o* V y y EFL 

y 

T 

s ® O' as‘s> + <o‘ aFS> - o ciss»^, 

O' 7 2 z t 

z 


C6 21> 


C6 22> 


C6 223 


Gomputei px’O^i’ammo usiri^ 2i>EGMENT 1 01 mulatiori Main t xjriGtlora 
of this pio^ramjrio i& to px'odiot tlma doperidfc?sit point coiiceutnatiori 
of a poilutaxit naiaa&od trom non-' & toady &oux ot; uridoi'" ti»ara.£.po3>t 
conditioris: Th« pr 0^1 ammo named "SEGMNT" ooriaiat ot various 

&ubpr ognarranes. 6 4> and its. lo^i-c is schiematically presented 

in Fig^ 6 5 The pnogi amme features ane similar to that of "PUFFS*’ 
fcfxoept tlie px ocedune for computing r ecepton concentration 
Further, if the i^eceptoi is not covered by the region of Influence 
<'SXR > XR> of the segmented plume or dispersion pax*an\eters 
evaluated at R^ point are negative or zex'O^ concentration 
contribution due to entire segmented plume is set to zero 

6 4 Mixed PUFF-SFGMENT for mulatlon 

The mixed FUFF- SEGMENT is a recently developed approach 
<:Zannettl, in which plume is descidbed by a series of 

elements^ may be segments or puffs^ or both The type of element 
affects computation of conceriti'ation field However^ element 

dyriamics remains unchang’ed CSec 6 3!!> This approach is more 

suited for ll^htwlnd conditions and can work in both calmwind and 
transport conditioris 

Identif iCtfitiori of element type The element type is a key 

t *Aotor in computiiiig' the plume coiiceritx''ation field The criteria 

for identifying the type of element is the ratio between its 

projected length on hr oizontal piano Ci e the horizontal 

distance between the first and second point of the eiement> and 

/y Pox a segment 
V 


D„/o > 2 

H y 


C6 24 J 
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Fig. 6.4 Progranme -tree for SEGMNT 








[STOP 


Fig. 6.5 


Schematic representation of the logic adopted 
in programme 'SEGMNT' 
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£ir\d^ for* a puff^ 

^ C6 25> 

H y 

Her*e <y v^ilue is. -taik^ri aus aivex'a^e of o* CIFS^ air&d CISS!) Since o* 

y y y y 

Gorifiriues fo wi*th fime^ all seg-merits wili event,ually become 

puff s. 

The efox'emenbioried cnitenie ais:s:ur*e& wht-u se^men-ts 

^ne -tr'ansfonmed into puffa^ the dlateinGe D between two 

H 

consecutive puffs will riot be eater than 20 * which is the 

V 

condition x*equined <Ludwl^ et o? ^ lyV7> tor -=i senies of puffs to 
provide an airriost perfect i epi eaentation of ^ coiitinuous. plirme In 
caimwind Gonditlons^ ^ 0 and puffs ai^e generated dinectly from 

tiic source 

Goncentx ation computation The concentr ation at a neceptor* R 
due to a cer tain source S most account fox* the contribution of all 
elements generated from S Thus the total conceritration at any 
I eceptor is the sum of the contnlbutlons of all the existing puffs 
plus the contribution of the closest segment 

A Puff Contribution* The concentration contidbution of a 
single puff at a receptor R is basically computed as indicated in 
Sec 6 2 However^ the c^s are evaluated at the centre point of 
the element CPuf f > as follows 

o* CIFS>+c/ <ISS> <y CIFS:>+o' CISS> 

P « y P 5 C6 2d> 

€/ 2 cr 2 

y z 

Pul f splitting is adopted whent-ver ^ The splitting 

procedure is same as in Sec 6 2 

B Segment Gontribution Contribution of only the closest 
se^ment^ which surrogates the effect ol all other se^ments^ is 
computed as described in Sec 6 3 

Computer pr o^x*amme for Mixed PUFF- SEGMENT formulation Main 
f unction of this programme is to predict tiino dependent point 
concentration of a pollutant released from non~steady source under 
li^htwlnd conditions The programme named “PUFSEG" consists of 
several subprogrammes <P1^ 6 6> and its io^ic is schematically 
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pieserit^ed in Fi^ 6 7 The bsts-lc str^uctur-e of the pr^o^r^imme is 
slmilaj> to thai of “'SE^MNT'* except foi^ the following chfiiri^es 

1 The plume consist of elements which may be segments or* 
puffs Thus it pequlT-ea identification of elements and seg-i'e^ation 
iri two diffeent senies of puffs and seg:ments 

2 It requires checking- coiitinuity ol segment 1 oi 

identif Ication of closest segment 

3 Though the pr ocedm e for a evaluation is same for* all 

elements Cpul f s or segments^ Fig6 7:>;^ the points of evaluation 
^re diffei’ent with i-ef enerice to physical location of the elements 
For puffs ar*e evaluated at centne point of the element while 

f or segments at R' location 

4 The concentnation contribution due to closest segment is 
computed first in subprogrammes SEOGON and to that is added the 
contribution of all the* puffs to get total point concentration at 
a receptor Segment coritributlon is computed only if the receptor 
comes in the influence region of the segment portion of the plume 
Cx coordinate of the flx*st point of the leading segment > XR> 

5 Tlie output gives complete plume profile specifying the 


type ol elements 
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r^r.nii 




li 


— *^ARE ALL RECEPTO 

►RS CONSIDERED ? 

Y 



■ARE METEOROLOGICAL CONDITIONS CHANGED ?j 


N 


lis STUDY DURATIOr^ OVER ? 


OUTPUT 'cTRIKlT /PLOT OF PLUME PROFICE“ 
RECEPTOR CONCENTRATION IN TABULAR FORM 


IStOP 


Fig, 6.7 


Schenatic represerTtation of the logic adopted 

In progranne 'PUFSEG'' 









7 PACKAGE OPERATION AND VALIDATION OF RESULTS 


7 1 Op^&r fiai^loriL 

The lo^diyi^ fiind compiling -Lo di&Ls. ot var^lous. pr-o^x^ammes: 
Lhixou^h TURBO PASCAL version 4 0 on Any IBM PC oompeLibie cr^eaLe 
files wiLh extension TPU/ EXE which can Lhen be executed in MS 
DOS envix onment wilhouL enbenin^ inLo TURBO Howevei^ sjome of Lhe 
TURBO PASCAL files <list ^iven in Appendix I> must be available on 

the disk dinectony of the package The execution of the package 
pr'Ognammes nequlnes a batch file which in itself can be considered 
as a master programme Ca set of DOS commands to execute various 
programmes constituting the package in a sequence detertulned by 
the user x*esporise]> A listing of INDISPOL batch file is presented 

in Appendix II The package execution stax’ts with a command 

OINDISPOL 

The package then displays first two cover pages CFlgs 7 1 and 
7 for convenlenoe all figures in this chapter are placed in a 

sequence at the end of the chapter > and a PAS^SGODE option CFlg 

7 3 ) The very first input the package requires is a package 

access code In the absence of the correct access code^ the 

package cannot be opened/initiallzed Once the correct code has 
been identified^ Inside cover* pages <lFlgs 7 4 and 7 5>^ general 

information about the package utilities CFlg 7 d a MENU HELP 

CFig 7 7> and SELECTION MENU CFlg 7 8> are displayed in a series 

The screen is divided in four windows The top most window gives 

heading^ the main window shows options available^ the side window 
displays selections made arid the bottom window indicates status of 
current operation arid options for forwardy^ackward movement within 
the package and transf erring control to the liard copy terminal for 
printing screen displays 
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The Relectlon of EDITOR option CPi^a. 7.8 And 7.0> 

^ovs Access to WORDSTAR 2000 <:WS2000> And dlsplAys EDIT OUIDE 
Exiting fpom the EDITOR mode T^etux'ns the iisei> to 
SELECTION MENU <P1^. 7.11>. 

User' is now needy to select Any slmulAtlon technlqpae 
<UiNIS^ MULTYS^ PUFFS SEOMNT on PUFSEF> depending upon the 
nequlnement. Howeven, the OUIDE option CF1^.7.12> enebies usen to 
know the feetunes of venlous simulAtlon tecliniques <Fl^s. 

7.13-7.17>. Aftei' the usen hiAS ^one thnoug^h the guidelines 
AVAiiAble^ peckA^e contnol is tnansfenned to the SELECTION MENU 
<F1^. 7.18> Again. Selection of a simulation technique displays 

its suboptions Ce.g. Fig. 7.19> until a tenminal selection 
<IFlg.7.20> is made. Usen can make multiple selections by choosing 
STAY <IF2!> option in which case the package netunns to SELECTION 
MENU. Multiple selections ane neconded in the sequence in which 
they ane selected in a dummy file which penmlts sequential 
execution of vanlous selections made. The selection NEXT CFig. 
7.20> nesults in execution of selected options. 


The package has pnovlslon of displaying stability 
classification table and coefficients adopted in computing 
dispex'sion standand devlatloris using powen laws CFlg.7.21>. 
Typical nesults of selecting UNIS- DIGITAL and UNIS-ORAPHICAL fon a 
specific input ane px'esented in Figs. 7.22-7.25. In a slmilan way 
othen simulation techniques cari be selected. The package pnovldes 
fon nepeatative selectioris aften all cunnent selections ane 
executed by tnartsfennlng contnol to SELECTION MENU^ if nequined. 
All digital outputs can be stoned in a usen specified files which 
can be accessed fnom within on outside the package. Exit fnom the 
package aften pnopen access nesults in display of end coven pages 
CFig.7.26 and 7.27!>. A sample Input/output fon vanlous simulation 
techniques adopted ane pnesented in Appendices III to V. 


7.2 Validation of Results 

Validation and llmltatioris of any simulation pnocedune 
muist be checked befone its application is made to a neal wox*ld 
pnoblem fon decision-making. The simulation pnocedunes wene 
validated tlinough a numben of tests. These tests Inclii^ed checking 


A 105908 

m • 1 i m brnm-n ,■ ^ - 


No. 
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stablll-ty of the sensltivl-ty of* the simulation pnocedunes 

to seveic £ii v£&nl£ables deflrdn^ the physios ot the pnoblem^ 
comparison of the results with other simulatiori procedures and 
compailson at the simulation results with experimentaiyf ieid 
measuj'ements In general predicted results were found to he quite 
logical and compared ver y well with those reported in earlier* 
studies CNema and Tare^ 1989^ Nema^ 1987> 




INDISPDL ; 

Interactive Package for Dispersion of Pollutants In the Atmosphere 
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F Ig, 7,2 Screen display of second ouf side cover page 


IND ISPD 
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FiQ. 7.3 Screen display of PASSCODE option 
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Fig. 7,6<a) Screen display of package general u-tlll"tles 
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<< CAPABILITIES OF UNIS >> 


This option IS built with the main objective to predict 
downwind, crosswind and vertical concentration profiles 
in stationary and homogeneous atmosphere using Gaussian 
model for a pollutant continuously released either from 
ground or an elevated point source at constant emission 
rate Input to this option can be through file "INPI’’ 
or terminal Digital output can be stored in a file to 
be specified when asked for 

Press <ENTER> to continue 

<< CAPABILITIES OF UNIS >> 


SAMPLE INPUT 

£ 00 i 1.00 1 00 

wind velocity(m/sec), at a height of(»), lapse rate (C/ 100m) 

1 

day <day*=1, means it is day? day=2, means it is night) 

2 

solar radiation is moderate (1 means strong, 3 means mild) 

1 6e-»-08 120 00 15 00 80 00 0 75 

source charact e r i s t i cs : 

mass emission rate (mi crogm/sec ), heat emission rate ( kj /se c . ) , gas exit 
velocity<m/sec ), dia of the eKit(m> 

0 50 

fractional reflection from the surface 

2 0 10 10 

Maximum, minimum and interval of downwind d i s tan ce s ( al 1 in km) 

Press <ENTER'> to continue 


Fig, 7,13 Screen display of guidelines for UNIS 
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<< CAPABILITIES OF MULTYS >> 


This option IS built with the main objective to predict 
pollutant concentral ions at selected locations ( called 
receptors) in a stationary and honriogene ous atmosphere using 
Gaussian model for a pollutant continuously released from 
ground or an elevated point source at constant emission 
rate Input to this option can be through file “INPP" 
or terminal Digital output can be stored in a file to 
be specified when asked for 

Sample INIPUT 

1 1 

No of receptor(s), no of source(s) 

Press <ENT£R> to continue 


<< CAPABILITIES OF MULTVS >> 


12.00 35 00 50 00 

receptor coordinates (km.m.m) 

2 00 11 .00 1 OO 

wind vel oc 1 ty (m/se c ) , height(m), lapse rate(C/100m) 

1 

day <day=1, means it is dayr day*2, means it is night) 

2 

solar radiation is moderate, (1 means strong, 3 means mild) 

0 00 0.00 O 00 

source coordinates <km,m,m) 

1 6e+08 120 00 15.00 60 00 0.75 

source characteristics* 

mass emassion rate<microgm/sec ), heat emission rate(kj/sec-), gas 
ve 1 o c i ty (m/s e c - ) , dia of the exit(m) 

0 50 

fractional reflection from the surface 

Press CENTER'^ to continue 


exit 


Fig 7.14 Screen ollsplay of gulolellnes for MULTYS 
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<< CAPABILITIES FOR PUFFS >>.... 


This option is built with the mam objective to predict time dependent 
point concentrations of pollutant released from fluctuating (non-steady), 
ground or elevated source at selected locations (receptors) in non- 
stationary and inhomogeneous conditions This option should be selected 
if wind velocity is low i e calmwind conditions (wind velocity should be 
<» 1 »/sec>- Some salient features of this option are as follows* 

1. Entire region of source influence is divided into meteorological grid 
and each block of the grid is allowed to take different meteorological 
parameters after every meteorological time step specified by you. 

2. Concent rat ion computations are carried out at the interval specified 
by you 

3. The plume profile specifies puff number, puff centre point coordinates, 
dispersion standard deviation and puff mass at a given time 

4 Delay time in getting the results directly depends on number of puffs 
generated, number of receptors and time interval at which con cent rat i on 
values are required So CAREFULLY SELECT these values 

b. INPUT to this option has to be through a file INP3 Output is in the 
digital form and can be stored in a file specified by you when asked for 


Hit <ENTER> to continue 


. « CAPABILITIES FOR PUFFS >> . 


However, you would soon get next version of the package which will allow you 
to draw i so concent rat i on lines 

6 Algorithm of this option is based on references Nema,A and Tare,V (1988) 
Computer Software for Atmospheric Dispersion of Pollutant”, presented at 
the ENVIROSOFT 88-2nd International Conference on Computer Techniques 
in Environmental Studies Porto Carras, Greece, September 19S8 and 
Zannetti, P (1988) ‘*New Mixed Segment-Puff Approach for Dispersion 
Modelling”, Atmospheric Environment, Vol 20, pp 1121-1130. 


Hit <ENTER> to continue 


Fig. 7.15Ca) Screen display of guidelines for PUFFS 
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<< CAPABILITIES FOR PUFFS >>. 

7 <SAMPLE 1NPUT> 

5^0000 ^0000 15000 30000 

Grid and block *xz©(«ll in m) XMAX , DELX , YMAX , DELY 
60 60 12 30 

Time var 1 abl ©s < a 1 1 in min ) 

Max time ,mel©orologi cal , element generation, and caculation time steps 
15 10 90 O 0 

Wind vel oc ity (m/sec ) , heighl(m), angle of velocity vectors from 2 and X axes 
1 50 0 0 0 0 200 

Lapse rate (C/100m) » precipitation rate(mm/hr), settling ve 1 o c i ty ( m /s e c ) , 

1 1 

Day<if 2 then it is night), strong solar r a d i at i on ( 2 , 3 meanmode r at e and mild 
1 

Dry de po 1 1 1 1 on ( 1 f 2 then it is wet deposition) 

Hit <'ENTER*> to continue 


<< CAPABILITIES FDR PUFFS >> 


0 0 1 

Scavenging ratio, rate constant ( /hr ) , minimum wind ve 1 o c i ty (m/sec ) 

15 120 

Gas exit veloci tyCm/sec ) , heat emission rate(kj/sec) 

1.6e-*‘08 I.Se-^OS 1.4e^08 1.3e+08 1 5e+08 
Mass emission rates ( mi c ro gm/sec) 

0.75 100 0.5 

Dia of exit(m), stack height(m) and fraction reflection from ground 
1 

Number of receptor(s) 

500 0 0 

Receptor coordinates (al 1 in m) 


Hit <ENTER> to continue 


Fig 7 15(b) 


Screen display of guidelines for PUFFS 
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<< CAPABILITIES OF SEGMNT >> 


This option is built with the main objective to predict lime dependent 
point concentrations of pollutant released from fluctuating (non-steady), 
ground or elevated source at selected locations (receptors) in non- 
stationary and inhomogeneous conditions This option should be selected 
for segmented plume i e transport conditions (wind velocity should be 
> H-3 m/sec). Some salient features of this option are as follows 

1 Entire region of source influence is divided into meteorological grid 
and each block of the grid is allowed to take different meteorological 
parameters after every meteorological time step specified by you 

2 Concentration computations are carried out at the interval specified 
by yo u 

3 The plume profile specifies puff number, puff centre point coordinates, 
dispersion standard deviation and puff mass at a given time 

4 Delay time in getting the results directly depends on number of puffs 
generated, number of receptors and time interval at which concentration 
values are required So CAREFULLY SELECT these values 

3 INPUT to this option has to be through a file 1NP3 Output is in the 
digital form and can be stored in a file specified by you when asked for 


Hit <ENTER> to continue 


<< CAPABILITIES OF SEGMNT >> 


However, you would soon get next version of the package which will allow you 
to draw 1 so concentrat ion lines 

6 Algorithm of this option is based on references Nema,A and Tart,V (1988) 
Computer Software for Atmospheric Dispersion of Pollutant", presented at 
the ENVIROSOFT 88-2nd International Conference on Computer Techniques 
in Environmental Studies Porto Carras, Greece, September 1988 and 
Zannetti, P (1986) "New Mixed Segment-Puff Approach for Dispersion 
Modelling", Atmospheric Environment, Vol 20, pp 1121-1130 


Hit <'ENTER> to continue 


Fig 7.16(q) Screen display of guidelines for SEGMNT 
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<< CAPABILITIES OF SEGMNT >> 


7. <SAMPLE INPUT> 

50000 50000 15000 30000 

Grid «nd block in m) XMAX , DELX , YHAX , DELY 

&0 60 12 30 

Tim« v«riable f ( al 1 in min.) 

Max t ime,mettorological , tlfiTient ganeralion, and caculation time steps 
1 5 10 90 O 0 

Wind ve lo c 1 ty (m/se c ) , he i ght <m ) , angl e of velocity vectors from 2 aind X axes 

1 50 0.0 0 0 200 

Lapse rate (C/100m) , precipitation rale ( frim/h r ) , settling ve 1 o c i ty < m/s e c > , 

2 1 

Night(if 1 then day)» sky is clear(if 2 means covered with cloud 
1 

Dry de pos 1 1 1 on ( 1 f 2 then wet deposition) 

Hit < ENTER > to continue 


<< CAPABILITIES OF SEGMNT >> 


0 0 1 

Scavenging ratio, rate constant(/hr), minimurri wind velocity<m/sec) 

15 120 

Gas exit ve 1 o c i ty < m/se c > , heat emission rate(kj/sec) 

1 6e+08 1 5e + 08 1.4e+08 1.3e-<“08 1.5e+08 
Mass emission rates(micro gm/sec) 

0 75 100 0 5 

Dia of exitCm), stack height(m) and fraction reflection from ground 

1 

Number of receptor(s) 

500 O O 

Receptor co o r d i nat e s ( al 1 in m) 


Hit <ENTER> to continue 


Fig. 7,16(lo) 


Screen display of guidelines forSEGMNT 
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<< CAPABILITIES OF PUFSEG >> .. 


This option It built with the main objective to predict time dependent 
point concentrations of pollutant released from fluctuating (non-steady ) , 
ground or elevated source at selected locations (receptors) in non- 
stationary and in-homogeneous conditions This option should be selected for 
moderate meteorological conditions ie light wind conditions (wind velocity 
« i-H m/sec) Some salient features of this option are as follows* 

1 Entire region of source influence is divided into meteorological grid 
and each block of the grid is allowed to take different meteorological 
parameters after every meteorological time step specified by you 

2 Concentration computations are carried out at the interval specified 
by you 

3 The plume profile specifies puff number, puff centre point coordinates, 
dispersion standard deviation and puff mass at a given time 

4 Delay time in getting the results directly depends on number of puffs 
generated, number of receptors and time interval at which concentration 
values are required So CAREFULLY SELECT these values 

5 INPUT to this option has to be through a file 1NP3 Output is in the 
digital form and can be stored in a file specified by you when asked for 


Hit <ENTER> to continue 


<< CAPABILITIES OF PUFSEG >> 


However, you would soon get next version of the package which will allow you 
to draw 1 s 0 c on cen t ra t i on lines 

6 Algorithm of this option is based on references Nema,A and Ta re , V . ( 1986) 
Computer Software for Atmospheric Dispersion of Pollutant", presented at 
the ENVIROSOFT 88-2nd Int ernat i onal Conference on Computer Techniques 
in Environmental Studies Porto Carras, Greece, September 1988 and 
Zannetti, P- (1986) "New Mixed Segment-Puff Approach for Dispersion 
Modelling", Atmospheric Environment, Uol.20, pp 1121-1130 


Hit <ENTER> to continue 


F,q. 7.17Ca) 


Screen display of guidelines for PUFSEG 
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<< CAPABILITIES OF PUFSEG >> 


7 <SAMPLE INPUT > 

5^0000 50000 1^^000 30000 

Grid and block fi 2 e(«ll in m) XMAX , DELX , YHAX „DELY 
60 60 12 30 

Time v*r 1 ables ( *1 1 in min ) 

Max 1 1 me , me teo r ol og 1 cal , e 1 tment generation, and caculation time steps 
1 5 10 90 O O 

Wind V e 1 o c 1 ty (m/se c ) , he 1 gh t (m ) , angl e of velocity vectors from 2 and X axes 
1 50 0 0 0.0 200 

Lapse rate (C/1O0m) , precipitation rate(mm/hr)r settling ve 1 o c i ty ( m/ s e c > , 

1 1 

Day(if 2 then night), strong solar radi at i on (2 , 3 mean moderate and mild) 

1 

Dry deposilion(2 means wet deposition) 

Hit <ENTER^ to ntinue 


. << CAPABILITIES OF PUFSEG >> 


0 0 1 

Scavenging ratio, rate c on s t an t ( / h r ) minimum wind ve 1 o c i ty ( m/se c ) 

15 120 

Gas exit vel oc ity (m/se c ) , heat emission rate(kj/sec) 

1.6e+08 1 5e-^0e 1 4e’^08 1 3e+08 1 5e-»-06 
Mass emission rates(niicro gm/sec) 

0.75 100 0 5 

Dia of exit(m)r stack height(m) and fraction reflection from ground 
1 

Number of receptor(s) 

500 0 0 

Receptor c oo rd i nat e s ( a 1 1 in m) 


Hit <ENTER> to continue 


Fig 7 17<b) 


Screen display of guidelines for PUfSEG - 2 
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Fig. 7.22 Screen display of digifal oufpuf In UNIS 
crosswind concentrai;ion profile 
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Fig. 7.23 Screen display of digital output in UNIS 
vertical concentration profile 
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Fig. 7 24 Screen display of graphical output In UNIS 
crosswind concentration profile 
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Fig. 7.25 Screen display of graphical outpul; in UNIS 
verilcal concentration profile 
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Fig. 7.26 Screen display of flrsf end cover page 
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8 SUMMARY 


Air pollu-tigirit^s which ar^e emiLled tx^am indiiiKhrlal and 
ohher man-mada sources aff'ec't hhe eco-system and its specific 
components^ leading to harmful consequences It^s^ therefox»ej, 
absolutely essential to stody the dispersion of pollutants for the 
envlrorimental assessment programme of polluted regions and to 
detei mine its impact on the environment For the estimation of 
pollutarits^ dispersion in the atmosphere several models are 
available To have simple and very uselul solutions of complex 
atmospheric phenomena^ computer aided slmuiatlon techniques are 
t^onsidex’ed to be very useful The Interactive Packages provide 
extx'a simplicity In Its use to the user In the present study 
efforts were made to develope a PC based Interactive Package to 
simulate pollutants^ dispersion from continuous or short-term^ 
steady or unsteady^ sm^face or elevated sources in stationary or 
non-statlonaj’y^ homogeneous or inhomogeneous atmosphere usin^ 
basic Gaussian theory The package corisists of a set of several 
programmes which are executed tfirou^h a batch control file in MS 
DOS environment usln^ TURBO PASCAL vex slon 4 0 The problems for 
dispersion of pollutarits are classified into five categories ^ UNIS 
foi sln^le^ steady^ and continuous source; MULTYS for steady and 
multiple continuous sources^ and PUFFS^ SEGMNT and PUFSEG for 
calmwlnd^ transport and ll^htwind conditions considering 
transformation and transpox tatlon of pollutants respectively 

Depending upon stack characteristics^ meteorological conditions 
and the physico-chemical nature of the eftiuent^ plume rise and 
hence effective stack height is calculated The power law is used 
for wind velocity profile in which power depends upon the 
lapse-r-ate Dispersion standard deviations ax^e determined by the 
mixiri^ or dispersive power of the turbulent flow witldn the 
planetary boundary layer UNIS provides the crosswind arid vertical 
profile of the pollutants assuming homogeneous and stationary 
atmosphere In the given domain In tabular and graphical forms 
MULTYS is exrtended form of UNIS in which there are multiple 
sources The principle of linear superposition is used to 



-the coiicent^pat^lon any nequined necap'tox Icoa-tioris. due 

i^Ci sevenal soux'ces ne^iect^in^ t,he int,e 2 >reneriGe be-tween any bwo 
PUFFS^ SEGMNT and PUFSEQ aimula'te pollubaribs^ dlspex'slon fx’om 
shon-t bei'm or eon-tlnuous^ sbeady oi l"luct,uat,in^ emissions ira non 
stal^ioriar’Y and inbomogeneous a-tmosphere The non-st^abionar-y 
Gondibions in mebeonoio^y and emissions ar-e bneabed as bhe poinb 

soui'oe emission inbo bhe abmosphene by a senies of insbanbaneous 

puf f &x^se^menbs/"eiemenbs whose char-acbenisbics ax^e updabed ab eaoh 
dispersion bime inberval Mebeonolo^ioai bhx^ee dimensional iTleldls: 
and emission papamebena ane also updLabed ab each mebeor-ologioal 

bime sbep For- bhe dynamics of a pufb/se^menb/elemenb^ boliowin^ 
bhiln^s axe corisldex^ed - ibs ^enenablon ab bhe sounce^ plume i*ise^ 
bi'ansponbablon by advecblve wlnd^ diffusion by abmosphenlc 


bur balance^ ground 

deposiblon which 

could 

be 

dx*y 

ox* 

web. 

arid 

ibs 

chemical 

bransf ormablon From 

bhe 

main 

MENU, which 

is 

specially 

designed 

fox bhls 

package^ bhe 

user 

can 


inbo 

any 

of 

bhe 

five 


menbioned cabe^onies of bhe problems The pticka^e provides on-line 
help> ^ulde and edibln^ facllibles which ai-e sell explanabony 

The package developed may pi'ove bo be a valuable bool 
1 ox esbimabin^ bhe exbenb of aln pollubion and decision making for* 
local and regional developirmenb However*^ bhei e is a need bo 

r*eflrie/modif y/e>iberid bhe package sbi*ucbur e/scope bo enhance lbs 

appllcabioris Following ane a few suggesblons on bhese lines 

1 The opbions wliich ane nob suppox*bed pnesenbly^^ bub 
ax'e available in bhe cunxenb INDISPOL sbr*ucbur*e^ should be fxxlly 
Implemenbed For* example^ bhe CONTOUR opbion can easily be 

developed by making pnovisloris for having access bo available 
package bo draw coriboux*s ^isoconcenbrablon llnes> on a similar* 
logic bhab is used for providing edlblng facillbles by giving 
access bo WORDSTAR 2000 

2 The compubablonal px*ocedux’e used in PUFFS SEOMNT and 
PUFSRO can be made more efflclenb by choosing varying puffy" 
segmenb/elemenb genex*abioriy"advecblon bime depending upon bhe 
calcuiabion bime and x ecepbox* locabions 



3 The s:lmiil£ihlori procedures t^heniselves should be modified 
<i> to ^ive statistioel computsitioris such as time aver^age 
conceritr'ation and standar^d deviatiori iii concentnatioris^ <lil-> to 
include effect of uppen layen inver^sion^^ concentr-ation computatior 
of secondary poilutarita^ vapiatiori in chemical i eactiori nate 
coristarit etc ^ arid <lii> to consider* abi upt changes in the 
meteoi'oiogical conditions 
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APPENDIX 1 LIST OF TURBO PASCAL FILES REQUIRED FOR 

INDISPOL EXECUTION 


GDRIVER TPU 
GKERNEL TPU 
GSHELLTPU 
GWINDOWTPU 


14-X9 FON 
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APPENDIX II INDISPOL BATCH FILE 


echo off 

path - cA.c:\rale«h;C:\v«2000. 

INTRODUC 

PERMIT 

If exist DUMY goto ter 
INFO 
toop 

ctol OPDUM 
del GUIDE • 
del DMLV * 

MAJNMENU 

if exist dumy goto ter 

RUNN 

If not exist ALC 11 goto IvH 
DOL1 

If exist DMLVII goto Ml 
If not exist DMLV1 goto arbitl 
EDIT 

c<J\wb2000 

W82 

cd cXrajeah 
arbitl 
CA03ABLE 
cIs 

UNI soon. 

.Ml 

If not exist ALC 12 goto Ivl2 
DOL2 


If exist DMLV12 goto M2 
If not exist DMLV2 goto arblt2 
EOT 

od\vwi2000 

wb2 

od c-Vajesb 
:Brbit2 
CALTABLE 
CiS 

UNISGRF 

:lvt2 

If not exist ALC. 21 goto Ivl3 
DOL3 

If exist DMLVId goto MO 
If not exist DMLV3 goto arbItS 
EDIT 

cxl\vwi2000 

¥A2 

cd c:\rajesh 
:arbit3 
CALTABLE 
cla 

RECPCONC 

:lv13 

If exist AUC.22 CX3NTOUR 
If not exist ALC.30 goto lvI4 
DOL4 

If exist DMLVH goto M4 
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If not CKi«t DMLV4 goto arWt4 
EDIT 

cd\v«2000 

va2 

od c-XraJeah 
:arb<t4 
CALTABLE 
cla 

PUFFS 

:lv!4 

If not «P(lat AlJC.40 goto Ivf6 
DOL6 

If «xlat DMLVIS goto M6 

If not epciat DMLV5 goto arbItS 
EDIT 

cd\wa2000 

v«2 

cd cNraJeah 
.erUtS 
CAOABLE 
cla 

SEQMNT 

:lv16 

If not CKlst ALC.60 goto Iv16 
DOL6 

If exiat DMLVie goto arbite 
If not exist DMLV6 goto arbltS 
EDIT 


cd\w8 2000 

W82 

od c:\rajeah 
:arbit6 
CALTABLE 
cia 

PUFSEQ 

:lvi6 

If exist QUIDE.61 QUID1 
If exist QUIDE.62 QU1D2 
If exist QUIDE.63 QUID3 
If exist QUIDE.64 QUID 4 
If exist QUIDE.66 QUIDS 
If exist QUIDE.* goto loop 
del ALC.* 

If not exist OPDUM goto notdon 
DOLL 

If not exist DUM1 goto notdon 
EDIT 

cd\\¥s2000 

wb2 

cd cXraJesh 
; notdon 

If exist OPDUM goto loop 
CHECK 

If exist DUM goto loop 

LAST 

:ter 
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APPENDIX III SAMPLE INPUT AND OUTPUT FILES OF UNIS 


Sfiimple InpurL File 


2 11 1 
1 
1 

1.7e+08 120 13 125 0.75 

0.5 

1.0 0.0 1.0 


Sample Oirtpcrt File 


concn. profile in y direction is 


X ( 1 n rri . ) 

Y < 1 n m ) 

Z ( 1 n 

1 000 

0 000 

156. 

1 000 

91 418 

156. 

1 000 

182.836 

156. 

1 .000 

274 254 

156. 

1 .ooo 

365 673 

156. 

1 000 

457.091 

156. 


m.) Concentration 

( mi c r ogm . / cum ) 


887 

138 826 

887 

126.61 1 

887 

96.044 

887 

60.600 

887 

31 .803 

887 

13.883 


concn profile in z direction is 


1 .000 

0.000 

156.887 

138.826 

1.000 

0.000 

125.510 

140.432 

1.000 

0.000 

94.132 

141.436 

1.000 

0.000 

62.755 

141.828 

1.000 

0.000 

31.377 

141 .603 

1 .000 

0.000 

0.000 

140.766 
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APPENDIX IV SAMPLE INPUT AND OUTPUT FILES OF MULTYS 

SAwipie Inpirt File 


10 

0 

50 

5 

0 

00 

0 

000 

0 

75 

0 

OO 

0 

000 

1 

OO 

0 

00 

0 

000 

1 

OO 

50 

00 

0 

000 

1 

OO 

-^0 

00 

0 

000 

1 

OO 

0 

00 

35 

000 

1 

OO 

0 

00 

75 

000 

1 

25 

0 

OO 

0 

000 

1 

50 

0 

00 

0 

ooo 

3 

00 

0 

00 

0 

000 

2 

50 

10 

OO 

1 

50 

1 

1 

0 

00 

0 

00 


0 00 


160000000 00 120 OO 15 00 1O0 00 0 75 

0 25 50 00 0 00 

100000000 OO 80 00 12 00 1O0 00 0 50 

0 25 -50 00 0 00 

100000000 00 85 00 20 00 75 00 0.50 

0 50 0 00 0 00 

100000000 OO 90 00 18 00 75 00 0 50 

1 00 O 00 0 00 

100000000 OO 40 00 15 00 40 00 0 50 

0 50 

SjBOYipie Ourtpixt File 


stability class* A which is most unstable 


no 

source coordinates 


eff ht discharge 


Kikm ,y*fri ,2fn 




m mi cr ogm /sec 

1 

0 00 0 00 0 00 



129 

56 1.6E-*-08 

2 

0 25 50 00 0 OO 



120 

71 1.0E-*-08 

3 

0 25-50 OO 0 OO 



101 

26 1 OE+06 

4 

0 50 0 00 O 00 



100 

61 I.OE-rOe 

5 

1 00 0 OO 0 00 



60 

09 1.0E+06 



pol lutant concn 

in 

desired receptors 


no receptor position 


c on cent ration 


KJkfTi ,y m , 

2 m 


micro gm /cu.m 


1 0.50 0 00 

0 

00 


612 874 


2 0.75 0 00 

0 

00 


1046.822 


3 1 00 0 00 

0 

00 


836 246 


4 1 00 50 00 

0 

00 


785 472 


5 1 .00 -50 00 

0 

00 


766 353 


6 1.00 0 00 

35 

00 


866 287 


7 1 .00 0 00 

75 

00 


870 685 


8 1.25 0 00 

0 

00 


1386.423 


9 1.50 0 00 

0 

00 


687 670 


10 3 OO 0 00 

0 

00 


27 866 
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APPENDIX V SAMPLE INPUT AND OUTPUT FILES OF 
PUFFS, SEGMNT AND PUFSEG 

Sample Input. Pile oP PUFFS, SEGMNT AND PUFSEG 


50000 50000 15000 30000 
60 60 2.5 30 
0.8 10 90 0.0 
1.50 0.0 0.0 200 
1 1 
1 


0 0 1 


15 120 


1 .6e+08 
1 .5e+08 
1 .6e+08 
1 .5e+08 


1 5e+08 
1 6e+08 
1 . 5e+08 
1 .6e+08 


1 .4e+08 
1 .5e+08 
1 .4e+08 
1 .5e+08 


1 .3e+08 
1 . 4e+08 
1 .3e+08 
1 .4e+08 


1 .5e+08 
1 .7e+08 
1 .5e+08 
1 .7e+08 


1 .2e+08 
1 .4e+08 
1 .2e+08 
1 .4e+08 


0.75 100 0.5 


10 


500 0 0 
1000 0 0 
1000 +50 0 
1000 -50 0 
1000 0 50 
1000 0 75 
1500 0 0 
2000 O 0 
3500 O 0 
4500 O 0 
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Sample Outpixt Pile oP PUFFS 


Plume profile *t 30 0 min. 


Puff 

No 

Centre point 
coo r di nates ( km ) 


Puff Mass 
( gm ) 

Rectangle 

s 1 gmay 
< (m) 

sigmaz ISC 

> Stb Cls 

1 

1.83 

0 00 

0 

19 

24000.00 

1 

1 

366 21 

1626.23 

1 

2 

1.67 

0.00 

0 

19 

22500 00 

1 

1 

336 30 

1330.27 

1 

3 

1 50 

0 00 

0 

19 

21000 00 

1 

1 

306 07 

1065 00 

1 

4 

1 33 

0 00 

0 

19 

19500 00 

1 

1 

275 48 

830 12 

1 

5 

1 17 

0 00 

0 

19 

22500 00 

1 

1 

244 48 

625 29 

1 

6 

1.00 

0 00 

0 

19 

18000 00 

1 

1 

213 01 

450 14 

1 

7 

0 83 

0 00 

0 

19 

22500.00 

1 

1 

180 97 

318.72 

1 

6 

0 67 

0 00 

0 

19 

24000 OO 

i 

1 

148 24 

209.94 

1 

9 

0 50 

0 00 

0 

19 

22500 00 

1 

1 

114 62 

124.08 

1 

10 

0 33 

0 00 

0 

19 

21000 00 

1 

1 

79 77 

61.53 

1 

1 1 

0 17 

0.00 

0 

19 

25500 00 

1 

1 

42 93 

22-88 

1 

12 

0.00 

0.00 

0 

19 

21000 00 

1 

1 

0 28 

29.23 

1 


Con cen t r at 1 on profile 

1 CRecp coordinates* 

at 30 0 
500,00 

min 

0 00 

0 

003 

c r ( mi c r ogm/ cu 

m) * 

475. 

, 13 

2 

CRecp 

coordinates 

1000 

00 

0.00 

0 

003 

cr (micro gm/cu 

m) X 

279 

70 

3 

CRecp 

coordinates * 

1000 

00 

50.CMD 

o 

003 

cr <mi c r ogm/cu 

flh ) * 

270 

76 

4 

CRecp 

coordinates 

1000 

00 

-50.00 

0 

003 

cr (mi c rogm/cu 

m ) X 

270 

76 

5 

CRecp 

coordinates * 

1000 

00 

0.00 

50 

003 

cr (mi cro gm/cu 

m > t 

285 

16 

6 

CRecp 

coordinates * 

1000 

00 

0.00 

75 

003 

cr (mi c rogm/cu 

m ) t 

286 

58 

7 

CRecp 

coordinates * 

1500 

OO 

0.00 

0 

003 

cr (mi c rogm/cu 

m ) X 

89, 

.63 

8 

CRecp 

coordinates 

2000 

00 

0.00 

0 

003 

cr (mi c rogm/cu 

m ) X 

24 

91 

9 

CRecp 

coordinates 

3500 

00 

0.00 

0 

003 

cr (mi c rogm/cu 

«rt) X 

0 

00 

10 

CRecp 

coordinates * 

4500 

00 

0.00 

0 

003 

cr (mi crogm/cu 

m) : 

0 

00 


Plume profile 60 0 min 

Puff Centre point Puff Mass Rectangle # s i gmay sigmaz ISC 

No coo rdi nates ( km ) (gm ) < (m) > Stb Cls 


1 

3 

83 

0.00 

0 

19 

24000.00 

1 

1 

708 

13 

7655. 

55 

1 

2 

3 

67 

0.00 

0 

19 

22500 00 

1 

1 

680 

54 

6974 

26 

1 

3 

3 

50 

0 00 

0. 

19 

21000 00 

1 

1 

652 

82 

6326 

03 

1 

4 

3 

33 

0.00 

0 

19 

19500 00 

1 

1 

624 

96 

5710 

70 

1 

5 

3 

17 

0 00 

0 

19 

22500 00 

1 

1 

596 

94 

5128 

14 

1 



97 


6 

3 

00 

0 

00 

0 

19 

18000 

00 

1 

1 

568 

76 

4578 

18 

1 

7 

2 

83 

0 

00 

0 

19 

22500 

00 

1 

1 

540 

44 

4060 

66 

1 

8 

2 

67 

0 

00 

0 

19 

24000 

00 

1 

1 

51 1 

93 

3575 

40 

1 

9 

2 

50 

0 

00 

0 

19 

22500 

00 

1 

1 

483 

23 

3122 

22 

1 

10 

2 

33 

0 

00 

0 

19 

21000 

00 

1 

1 

454 

32 

2700 

93 

1 

1 1 

2 

17 

0 

00 

0 

19 
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